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[ Abstract. 

D ■ The peculiar effects of a quantum measurement are completely foreign to classical 

f"*^ , physics, but relevant features of an entangled quantum state are accessible to 

semiclassical analysis. Classical surfaces, or manifolds in phase space correspond to 
quantum states in Hilbert space. Subsystems specifie factor spaces of the Hilbert 
space, which correspond to lower dimensional phase spaces. An entangled state 
corresponds semiclassically to a surface that cannot be decomposed into a product 
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. of lower dimensional surfaces, for a specific choice of subspaces. Such a classical 

' factorization never exists for ergodic eigenstates of a chaotic Hamiltonian. 

Entanglement is best described in the language of density operators. The product of 
\q , the bra space and the ket space forms a dyadic basis for quantum operators, analogous 

to combining the Hilbert spaces for subsystems. Thus, the space of quantum operators 
. corresponds to a double phase space. The various representations of the density 

\\ operator then result from alternative choices of allowed coordinate planes in this 

enlarged phase space. The chosen plane may be a phase space on its own, as in the 
case of the Wigner function, akin to a probability distribution. Its Fourier transform, 
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^ ' the chord function, or the quantum characteristic function, lies in an alternative phase 

O^ 1 space. 

The reduced Wigner function represents the partial trace of a density operator 
over a subsystem and so contains contains all information about a given subsystem. 
It results from the projection of the original Wigner function, just as in the definition 
of a marginal probability density. The reduced chord function is a section of the 
original chord function and also contains all relevant information. The purity of the 
reduced density operator, i.e. the square of its trace, is a good measure of entanglement, 
obtained by integrating either the square of the reduced Wigner function, or the square- 
modulus of the reduced chord function. 

The Wigner function can be associated to probabilities for measurements of the 
eigenvalues, ±1, of the operators corresponding to reflections in phase space, i.e. 
generalized parity operators. These are wavelike properties which appear unintuitive if 
the quantum system is pictured as a particle. Bell inequalities for such measurements 
can be violated even for classical looking states with positive Wigner functions that 
have evolved classically from product states. These include the original EPR state. 

Entanglement with an unknown environment results in decoherence. A simple 
example is that of the centre of mass of a large number of independent particles, 
entangled with internal variables. In this case, the Central Limit Theorem for Wigner 
functions leads to some aspects of Markovian evolution for the reduced system. 
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1. Introduction 

The realization that quantum mechanics admits entangled states goes back to 
Schrodinger in 1926 pQ. He went on to coin the term entanglement in 1935 [2], but 
the dramatic example of a biological cat coupled to a decaying nucleus was never meant 
to be operational. Einstein, Rosen and Podolsky (EPR [3]) discussed an example of a 
simple entangled bipartite state in the same year. Their concern was the compatibility 
between Heisenberg's indeterminacy principle and the generation of strong correlations 
through a measurement on a member of a pair of particles, even when they could no 
longer be interacting. It was the formulation of Bell inequalities, starting in 1964 [I] 
that provided a litmus test for nonlocal correlations in quantum mechanics. The initial 
concern was centred on hidden variable theories and the possibility of their emmulating 
quantum correlations even for particles that have ceased to interact. Such violations 
of local causality, detected by Bell inequalities, could not have developed within any 
kind of classically evolved ensemble, irrespective of whether the variables are explicit or 
hidden. 

Quantum information theory [6] has given a new boom to the study of the 
qualitative distinctions between classical and quantum mechanics and to establishing 
their quantitative measures. There is nothing so dramatic about the development of 
nonclassical correlations between particles that are still undergoing an interaction, but 
this question has acquired promissing applications in future quantum computations. 
Necessarily, these deal with finite dimensional (Hilbert) state spaces, for which the 
appropriate entanglement measures are now well established. 

One of the difficulties in aplying semiclassical methods to the study of entanglement 
is that the former have been developed for infinite dimensional Hilbert spaces. Not only 
are these an extrapolation from the few qubits that have been usually considered in 
quantum information theory, but entanglement is most clearly exhibited through the 
correlations in elementary either- or experiments. This seems to privledge simple state 
spaces of a single qubit, such as spin-1/2 systems. For this reason, Bohm's version of 
EPR |7J has become much more popular than the original full phase space version. A way 
around this difficulty is to consider the measurement of special observables which have 
only a pair of eigenvalues, even though they operate on states within an infinite space. 
It turns out that one of the most renowned phase space representations in quantum 
mechanics, the Wigner-Weyl representation, is based on such operators. Usually, this 
representation is viewed as a way of eliciting classical features in a quantum state, but 
it will be used here mainly as a probe into nonclassical correlations. 

The development of semiclassical theory throughout the last century allows us to 
trace the classical skeleton underlying many features of quantum evolution. These 
classical structures are the core of approximations that improve asymptotically in 
the limit of large classical actions, or, more formally, as Planck's constant, h — > 0. 
In the case of a finite dimensional Hilbert space, this becomes the limit of large 
dimensions. Even though entanglement is a subtle phenomenon, it leads to gross 
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Figure 1. The classical trajectories, xi(t) and X2(t), for each component in its own 
phase space, are projections of the full trajectory for the entire system. 

violation of inequalities and to quantitive measures which are not beyond the accuracy 
of semiclassical approximations. Therefore, it is appropriate to enquire into the manner 
in which classical structures can be implicated in such a very nonclassical feature of 
quantum mechanics. 

Traditionally, semiclassical theory is concerned with the unitary quantum evolution 
of closed systems, which are thus described classically by Hamiltonian dynamical 
systems. Each point in phase space accounts completely for the state of the classical 
system which evolves along a trajectory. A bipartite, or multipartite system is 
accomodated in this correspondence by a higher dimensional phase space. Each point 
still evolves as a single- dimensional (1-D) trajectory, but its projections onto the 
subspaces, which describe the succesion of possible states of each component of the 
system, are also 1-D trajectories in their own right, as shown in Fig. HJ 

Part of the power of Hamiltonian dynamics lies in the freedom to transform between 
different sets of phase space coordinates. This canonical invariance emphasises the 
importance of the unified evolution of the full system over that of the component 
trajectories, which are seen to depend on the particular choice of coordinates. In 
contrast, it is in the separation into components that the phenomenon of quantum 
entanglement emmerges. The particular nature of quantum measurement lies behind 
this difference, as is discussed in section 2: It is only when this is combined to the 
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preceding unitary evolution, that the unclassical correlations between the components 
become manifest. 

Therefore, the study of the properties related to entanglement should be viewed as 
an objective that is imposed externally on semiclassical physics, which perhaps explains 
the low priority received by this goal so far. In these lectures, we will only be concerned 
with the most elementary kind of entanglement, i. e. that of pure bipartite states, for 
which the measures of entanglement are well established. Even so, it will be seen that 
this simple case requires the introduction of theoretical instruments of semiclassical 
theory that are far from elementary. Not only do we need to cope with a higher 
dimensional phase space for the description of a biparite system, but it will be shown 
how the simplest semiclassical description of operators is achieved in a phase space 
with double the dimension of the one corresponding to the states on which they act. 
Conversely, some of the most relevant structures for entanglement, such as partial traces 
and probability densities, can be interpreted as projections of the Wigner function, or 
sections of its Fourier transform. 

The following section reviews the different ways in which features of quantum 
mechanics, interference and entanglement, are nonclassical. Simple examples introduce 
the reflection symmetries, quite familiar for classical waves, that will play a major role 
in the Wigner function formalism. Prior to this though, it is useful to consider classical- 
quantum correspondence in a more simple-minded way. This is the subject of section 
3, which introduces product spaces for both quantum states and classical probability 
distributions in phase space. In either case, the factorizability is broken by an interaction 
Hamiltonian, leading to correlations for measurements on the different components. In 
the classical case, these correlations are constrained by general Bell inequalities. This 
section also introduces the Schmidt decomposition of quantum states. 

Section 4 reviews standard semiclassical theory for quantum states. Special 
emphasis is given to products and factorization of both the phase spaces themselves 
and the internal Lagrangian surfaces that support the quantum states. This product 
structure is then generalized in section 5 to the representation of operators. Dyadic 
operators of posistion or momentum eigenstates form a complete linear basis for all 
quantum operators, which correspond to planes in double phase space. Linear canonical 
transformations take these into the phase space coordinates for the Weyl representation 
and the chord representation, its Fourier transform. These bases are associated 
respectively to phase space reflections and translations and to the corresponding 
quantum operators. In the case of the density operator, we thus obtain the Wigner 
function and the chord function, both presented in section 6. Though the Wigner 
function cannot be interpreted as a probability distribution, because it may be negative, 
it coincides with the difference for probabilities of measuring either the positive, or the 
negative eigenvalue of the reflection operator. Section 7 is dedicated to projections of 
the Wigner function and sections of the chord function, which represent the reduced 
density operators. The loss of purity of the latter, obtained as integrals of the square 
of either the reduced Wigner function, or the reduced chord function, indicate that the 
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overall state is entangled. 

It may be guessed that an initially classical pure state, the product of Gaussian 
Wigner functions, would not be entangled by a simple rotation of positions and 
momenta. After all, this class of states, including the original EPR states, could 
stand in for a classical phase space distribution. However, this is not so, as shown 
in section 8: The reflection correlations for such states violate Bell inequalities, even 
though measurements of positions and momenta can only correlate classically. The 
transformation to centre of mass coordinates for any number of particles, studied in 
section 9, has similar features. By invoking the Central Limit Theorem for Wigner 
functions, we obtain features of the nonunitary evolution of the centre of mass in 
agreement with Markovian theory, i. e. the exact solution of the Lindblad equation 
for the density operator. 

The final section relates double phase space geometry to the semiclassical Wigner 
and chord functions. These are not known in detail for eigenstates of chaotic 
Hamiltonians, but it has been proved that ergodic eigenstates are supported by the 
entire energy shell. In this case the unitary transformation which factorises the state 
can have no classical correspondence. 

A lot of the experimental work related to entanglement has been carried out in 
quantum optics. Rarely is the full generality of semiclassical states employed there and 
one can rely mainly on states derived from the eigenstates of the harmonic oscillator, 
even when phase space is invoked [8]. For this reason, the initial examples of phase 
space structures are here chosen among states of this type. The reader who wishes to 
avoid the more subtle aspects of semiclassical theory can mostly skip sections 4, 10 and 
parts of 5. 

2. Entanglement and classical physics 

Entanglement is considered to be a quintessential quantum property which defies all 
attempts at a classical correspondence. For this reason, its description in terms of 
the classical concept of phase space might appear foolhardy. It could be that the 
semiclassical program of uncovering meaningful relationships between XlX'th and XX'th 
century mechanics would be overstretched. Perhaps, though, such an endeavor would 
make more sense if it were recalled that the usual validity of a classical description of 
macroscopic phenomena can be attributed to the effect of decoherence. In its turn, 
this results from the entanglement of a given system with an unco ntro led environment, 
caused by interactions that can be minimized, but never entirely eliminated. Thus, in 
spite of the fact that the common working languages employed in classical and quantum 
mechanics are quite alien to each other, it is hard to fully comprehend why the outcome 
of decoherence should be the emerging appropriateness of a classical description for 
quantum systems, unless we can detect its traces even within entanglement itself. A 
simplified version of this program will be sketched in section 10. 

Before attempting to establish a bridge between some features of quantum 
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entanglement and classical mechanics, it is worthwhile to consider the more obvious 
way in which interference already separates these theories. In contrast, the analogy 
of quantum mechanics with classical waves is much smoother: The latter may be 
superposed linearly and they interfere in the same way as matter waves. In a simple 
two-slit experiment, the initial quantum state is prepared as a coherent superposition 
of momentum eigenstates, with eigenvalues that can be classically measured: The 
probability for each momentum direction is the same as for a uniform ensemble of 
classical states. The evolution through a pair of slits generates classical interference, 
equally observable in water waves, or sound waves. Quantum strangeness only emerges if 
the intensity of the resulting interference pattern for the conjugate variable, the position, 
is interpreted as the probability for the position measurement of a single particle, 
moving according to classical mechanics. Even so, the particular nature of quantum 
measurement itself does not play a prominent role in the phenomenon of quantum 
interference. The subsequent quantum state is certainly redefined by the result of the 
measurement, but this is not a crucial feature of quantum interference, no matter how 
unclassical its interpretation for a single particle. 

The success of the semiclassical treatment of interference phenomena is no real 
surprise. If we start from Feynman's path integral formalism [TOlfTTj . quantum evolution 
is described by a continuum of interfering paths. Semiclassical theory merely groups 
these around a few particular classical trajectories with their Feynman phase. The 
amplitude of each of these discrete interfering terms is then given by a local integration 
over the continuum of paths. Classical mechanics takes its part in the theory, but there 
is no limitation to classical phenomena and interference is well described. Indeed, the 
role of classical mechanics is the same as ray optics in classical wave theory. 

In contrast to interference, the unintuitive nature of entanglement is derived from 
that of quantum measurement itself. In no way does this tally with the common sense 
description of the macroscopic world. Nothing in our everyday experience prepares us 
for the collapse of a state that is measured into one of several possibilities. The common 
sense presuposition would be that the effect of the measurement on the system could 
and should be made negligible. Entanglement highlights this phenomenon in a specially 
subtle way, because it involves pairs of measurements on systems with at least two 
degrees of freedom, or components. 

If we consider classical waves, or particles, it would be indeed strange to imagine 
that such a collapse could result from the measurement of a subsystem, thus constraining 
the possible states of the complementary subsystem: It is well known that playing a 
note on a piano, i.e. exciting a finite string, will provoke a response on the next octave 
string. Here we have two nearly independent systems, stretched strings, weakly coupled 
by the surrounding air. Perhaps, it is better to consider the same note on two nearby 
pianos, so that we consider the interaction of identical systems. The wave form assumed 
instantaneously by the pair of strings may be used to describe the state of the whole 
system, or else, we may prefer the Fourier representation, in terms of the eigenstates 
for the discrete set of allowed frequencies of each string. These classical strings are 
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completely analogous to the textbook example in quantum mechanics of particles moving 
in 1-D, each in its own box. But there is no way in which a photograph of one of the piano 
strings will affect the sound produced by the other string, no matter how entangled the 
quantum analogues happen to be! Likewise, the measurement of the frequency spectrum 
of the vibrations of one of the strings does not oblige it to choose among the various 
overtones and we would be even more surprised if this led to a correlated jump in the 
other string. 

Yet this is just what we would expect for the analogous quantum system composed 
of two particles in their 1-D boxes, coupled by the same Hamiltonian that may account 
for the atmospheric interaction. Such a measurement would single out a discrete energy, 
or equivalently a discrete momentum modulus. Furthermore, in the quantum system, 
we could also measure the position of the particle, with a probability density that is 
specified by the wave intensity. No equivalent interpretation can be imputed to the 
classical wave, so that such a position measurement would then be devoid of meaning. 

@ 

Just as there are measurements on a quantum system that are meaningless for a 
classical wave, there are others which make no sense for a classical particle. Consider 
the excitation of a piano string by the same note, but played on a clarinet. This has 
only even harmonics, because it is equivalent to a string that is free on one side. Then 
only the even harmonics will be excited in the string, which will be symmetric about its 
midpoint. Such an even (or odd) parity, i. e. the symmetry (or antisymmetry) of the 
classical stationary wave is certainly a measurable property of the analogous quantum 
state. Indeed, even a classical wave, a string that is free on one side, could in principle be 
used as a probe to measure directly the even component of the wave, instead of exciting 
it. But what would it mean to measure the parity of the corresponding classical particle 
in a box? Generalizations of such parity measurements, distinguishing the eigenvalues 
of non-mechanical observables, will play a major role in the following discussions of 
entanglement. 

Measurement theory lies outside the scope of a semiclassical treatment. However, 
such experimental outcomes will be preceded by (unitary) quantum evolution, which 
is not so adverse to a classical description. Indeed the process by which subsystems 
become entangled is a preparation that precedes any quantum measurement. It is only 
in the probabilistic interpretation of the subsequent measurement on the system that the 
quantum and the classical viewpoints fundamentally diverge. As it happens, standard 
measures of entanglement require that the components of the system be completely 
defined, but do not pre-specify the measurements to be performed. Thus, the presence 
of entanglement only indicates the possibility that some subset of measurements will 
have nonclassical correlations. It is precisely this lack of definition with respect to future 
quantum measurements that allows space for a semiclassical treatment. 

The study of classical waves displays many of the properties of a simple quantum 

§ It should be remembered that the classical particle analogy here is not related to the phonons that 
are generated by second quantization within each mode. 
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system. Indeed, Rayleigh's The Theory of Sound [12] anticipates some results later 
rediscovered in semiclassical theory. However, each piano string is a system with infinite 
degrees of freedom. Though it is not forbidden to consider coupled fields [jj, the following 
lectures will concern mainly systems with a finite number of degrees of freedom. In most 
cases, two degrees of freedom already suffice to discuss the relation between the concept 
of entanglement and classical mechanics. So we start with a review of classical-quantum 
correspondence. 



3. Classical-quantum correspondence 

The simplest quantum systems with a classical correspondence have a single degree of 
freedom, e.g. a particle constrained to move in a straight line. The classical state of the 
system is described by its position, q, and its momentum, p. Together they define a point 
in phase space, x = (p,q), which is a 2-D plane. Perhaps, classical state space would 
be a more appropriate term, because each point specifies all future motion of a classical 
system, once the Hamiltonian, H(x), is specified, through Hamilton's equations: 

P=-— q = — (3 1) 

dq ' dp 

These equations may be compactified into the form 

± = 3§ , (3.2) 
with the definition of the (2 x 2)-dimensional matrix 

(3.3) 



" 


-1 " 
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acting on the phase space points, x = (p,q). Unless H(x) is quadratic, this motion is 
nonlinear. 

Corresponding to this 2-D plane, quantum mechanics matches the states \ip) of 
an infinite dimensional Hilbert space, H, on which act the operators, q and p. Each 
eigenstate of q, labled by the eigenvalue q , corresponds to the vertical line, q = q , 
whereas the horizontal phase space lines are matched by eigenstates of p. These 
operators do not commute, [p, q] = iH, but if we appropriately symmetrize the order 
in which p and q appear in H(x), then the motion of the states \ip) is also determined 
by the quantum Hamiltonian H(x), through the linear equation, 

ih^W) = H{x)\1>), (3.4) 

i.e. Schrodinger's equation. 

The uncertainty principle excludes the existence of a quantum state that 
corresponds precisely to a phase space point. However, the unavoidable dispersion 
in measurements of position, or momentum allow us to seek an approximate 
correspondence with probability distributions of phase space points. This is 

|| Perhaps, quantum superstring theory will tackle entanglement someday. 
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unsatisfactory as far as interpretation is concerned, because probabilities are matched 
to the square of a state rather than the state itself. Nonetheless, a certain intuition can 
be obtained through this analogy. Given a phase space probability density, f(x), the 
expectation value of any classical observable 0(x) is given by 

E(0) = J dxO(x) f(x). (3.5) 

Hence, the dispersions in position and momentum are 5q 2 = E((q — E(q)) 2 ) and 
Sp 2 = E((p — E(p)) 2 ). The uncertainty principle then imposes that only phase space 
distributions for which A' = SqSp > h should be considered. 

However, this quantity is not a classical invariant. The flow, x(0) — > x(t), generated 
by the Hamiltonian is a canonical transformation, so that |13| . for all t, 



p(0)-dg(0)= & p(t)-dq(t), (3.6) 

'70 Jjt 

where 70 is any circuit and 70 — > 7*. General Hamiltonian evolution will stretch and bend 
any closed curve that is initially compact, so that a probability distribution that is unity 
inside 70 and zero outside will not have constant A'. Linear canonical transformations, 
that is, symplectic transformations, are well known to be specially favourable for 
classical-quantum correspondence, as will be further discussed. It will be shown in 
section 6 that, A, the determinant of the covariance matrix, 



K 



5p 


5pq 


5pq 


5q \ 



(3.7) 



where {bpq) 2 = E(pq — E(p)E(q)), is invariant under symplectic transformations. 

To discuss entanglement, we need more than one degree of freedom. Quantum 
states can then be decomposed into a basis of product states, 

MHV'i) «>••#£>, (3-8) 

which span the full Hilbert space, H = Hi <%> ...Hi <%> ...Hl, i.e. the tensor product of 
the factor Hilbert spaces that describe each degree of freedom. Likewise, the full phase 
space is now a Cartesian product of 2-D conjugate planes, each the phase space for a 
particular degree of freedom, 

x = x\ x ...x\ x ...xl (3.9) 

and thus has 2L dimensions. However, we must be wary of the difference between the 
classical and quantum geometries: a phase space strip, 6q, corresponds to this range 
of eigenvalues for the operator q. This set of eigenstates spans an infinite dimensional 
subspace of the product Hilbert space, whatever the number of degrees of freedom. On 
the other hand, each of these position eigenstates corresponds to one of the parallel L-D 
g-planes within the 2L-D phase space strip. 

The classical or quantum motion for systems with more than one degree of freedom 
is still defined by a Hamiltonian, H(x), or H(x), but now dH/dx is a 2L-dimensional 
vector and J is a block matrix. We shall also use the skew product, 

L 

x A x' = ^{piq[ - qip'i) = J x ■ x' . (3.10) 

n=l 
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This symplectic area of the parallelogram formed by the vectors x and x' is 
invariant with respect to symplectic transformations. Again, these are linear canonical 
transformations, with ( 13. 6p interpreted as a line integral in the 2L-D phase space. For 
higher dimensional systems, all even dimensional volumes, from 2 to 2L, are preserved 
by canonical transformations [T3] . 

If the degrees of freedom are completely decoupled, each with its own probability 
distribution, fi(xi), the full probability distribution will be just the product, 

In this case, the probability distribution for a single degree of freedom is reobtained by 
tracing over the other variables: 

fi{xi) = / f{x) dx 2 ...dx L . (3.12) 

If the full probability distribution cannot be factored into a product, the above equation 
then defines the marginal distribution. This process foreshadows that of partial 
tracing over the density operator, to be studied in section 7, which is central to the 
study of entanglement. 

The product nature of the classical distribution will be retained throughout the 
evolution, if the full classical Hamiltonian is purely additive, 

H(x) = H x {x x ) + ..Jf,(z,) + ...H L (x L ), (3.13) 

i.e. if there is no coupling between the motions of the several degrees of freedom. This 
follows from the decoupling of Hamilton's equations into 

±i = J,^i (3.14) 

OXi 

for each degree of freedom. In other words, if k ^ I, then xi(t; xi ) does not depend on 
Xk (nor on the initial value, x^q). Furthermore, we then have 

fi(t; xi) = fit; xi{-t; Xi)) and f(t; x) = fi(t; Xi)...f L {t; x L ), (3.15) 

where Xi(—t; xi) specifies the past location of X[. Likewise, the volumes in each subspace 
will be preserved and the conservation of the 2L dimensinal volume is just that resulting 
from the conservation of the factor volumes. 

For a classical system, the transition from product probabilities to general 
probabilities can only be generated by coupling terms in the driving Hamiltonian, 
containing cross products, which are at least bilinear in the different variables. A 
general classical observable will be a function of all the phase space variables and its 
expectation is accordingly given by (I3.5p . For instance, this might be the either-or 
observable, 0\ = ±1 for detecting some physical property associated with one particle, 
or the detection of 2 = ±1 for a second particle. For classical particles which have 
been allowed to drift sufficiently far from each other after interacting, the result of the 
Oi-measurement will not affect the (^-measurement and vice versa. Therefore, the 
correlation must be represented in the form 

£7(0i ; 2 ) = £(Oi O z ) = J d Xl dx 2 Oi(xi) 2 {x 2 ) f{x u x 2 ). (3.16) 
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This equation has the same form as correlations postulated for local hidden variable 
theories [H [H] . Indeed, one of the reasons for this choice is that (13.161) must hold for 
any evolution of f(x) governed by classical mechanics. This form for the correlation 
between different components of the system is then taken as a prerequisite for theories 
that in all other respects should give the same results as quantum mechanics. Since 
this is certainly not one of the objectives of classical mechanics, such conjectures then 
necessarily demand extra, unknown and hence hidden variables. 

It is due to the seminal work of Bell pE] that we are able to compare, through 
inequalities, the correlations predicted by quantum mechanics with a very wide range of 
possible local correlations. The point is that any measurement affects the entire quantum 
state, i. e. both its components, unless the state happens to be an eigenstate of the 
measured observable. So quantum measurements are not local in the sense that led to 
( 13.16!) . In case of the general CHSH inequality pj2,[6], involving either-or observables, 
0\a, Oib, 2a and 2 f>, (13.161) implies that 

\E(0 la ; 2a ) + E(0 la , 2h ) + E(0 lb ; 2a ) - E(0 lb ; 2b )\ < 2. (3.17) 

As well as constraining possible hidden variable theories, this inequality can be used as a 
detector of nonclassical correlations in quantum mechanics. This kind of nonclassicality, 
entanglement, is much more subtle than quantum interference effects, as will be 
discussed in the later sections. A dip into Bertlmann's socks and the nature of reality 
[9] by Bell provides a delightful discussion of all the main points concerning classical 
locality versus quantum correlations. The book by Peres [H] is also recommended. 

It is worthwhile to discuss some specific examples of systems with more than one 
degree of freedom. An obvious possibility is a collection of particles, each moving in one 
dimension. Another is a single particle moving in two, or three dimensions. Classical and 
quantum mechanics make no distinction between these alternative interpretations of the 
dynamical variables. All that is demanded is that the variables pertaining to different 
degrees of freedom commute, \pk,Qj] = ifiSkj, or > correspondingly, that the classical 
Poisson bracket {pk,qj} = 5kj (see e.g. [16]). We can also use angular momentum and 
their conjugate angles. But are other variables, obtained through classical canonical 
transformations, allowed? 

For example, consider our piano string, now modeled as L masses connected by 
harmonic springs. We can switch to the L normal modes of vibration. This is a linear 
canonical transformation, which substitutes the original L conjugate planes, x\ = (pi, qi), 
by new conjugate planes, x\ = {jp'i,q[), that now describe collective motions of the L 
masses. This is also a proper phase space to be quantized, x\ ^x\. Another important 
example of a quantizable canonical transformation follows from the description of a 
collection of particles in terms of the centre of mass together with internal coordinates. 

Whatever the physical realization, symplectic transformations, correspond exactly 
to unitary quantum transformations and hence to equivalent quantum systems [17] . 
These transformations generally redefine the components of the full system and may take 
an entangled state into a product state, or vice versa. Any measure of entanglement 
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is affected by such a general transformation, so one requires only that the measure 
be invariant with respect to local unitary transformations, lying within each separate 
component. As for nonlinear canonical transformations, these are not exactly matched 
by quantum unitary transformations [IT] and, hence, cannot be directly quantized. It 
might still be useful sometimes to push this correspondence through, but it must be 
remembered that the result is only a semiclassical approximation. 

Taking again the continuum limit, L — > oo, each normal mode of the finite chain 
converges onto one of the lower modes of the continuous string. There is no essential 
difference between the interaction and hence the entanglement among these modes of 
the continuum and that of finite modes (caused by residual nonquadratic terms in the 
Hamiltonian). In each case, there corresponds a plane in the phase space, which is 
of infinite dimension in the case of a field. The entanglement between modes of the 
electromagnetic field within a finite cavity also has a similar interpretation in terms of a 
classical field. The unperturbed motion is now that of a quantized harmonic oscillator, 
corresponding to a classical oscilation in each phase plane . 

Another example is that of a particle with internal structure. The latter may be 
described by an angular momentum, coupled to the translational degrees of freedom 
by an external field. The Stern-Gerlach experiment is just such a system, in which the 
magnetic moment, tied to the spin angular momentum of the electron is coupled to 
its position by an inhomogeneous magnetic field. The spin is an intrinsically quantum 
mechanical two level system and the interest in quantum information theory tends 
to emphasise such simple quantum systems. But, in principle, there is no difference 
between this case and a Rydberg atom, prepared in a state with a large electric dipole 
moment, coupled to position through an inhomogeneous electric field. Such a system 
can be described more naturally in classical terms. Cavity quantum optics deals with 
the coupling and hence the entanglement of the internal states of individual Rydbreg 
atoms with a specific mode of the electromagnetic field. 

For all these systems, coupling terms in the overall Hamiltonian will destroy the 
product form of an initially decoupled quantum state, or classical distribution. We 
should bear in mind three basic differences between classical and quantum systems: 
i) the nature of the initial state; ii) the nature of the evolution and iii) the effect of 
experiments. As we have seen, the last is the most radical difference, which, indeed, 
gives rise to the concept of entanglement. Our objective here is to cast the quantum 
mechanical description of i) and ii) in the most classical terms possible, so as to highlight 
the truly inovitative elements of the quantum theory when iii) is considered. 

A fundamental difference between the quantum and classical descriptions should be 
discussed before proceeding: The analogy between the evolution of classical probability 
distributions and quantum states is somewhat deceptive in as much as the latter 
determine only probability amplitudes which can be complex and interfere with each 
other. To arrive at a closer analog of probabilities, we should, in some sense, square 
the quantum states. The correct procedure is to define density operators, or their 
phase space representation, Wigner functions, to be studied in section 6. However, their 
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evolution is nonclassical, unless the Hamiltonian is quadratic. 

It will be only in the context of the density operator that it becomes meaningful 
to distinguish between pure states and mixed states. Taking an average over a set 
of probability distributions defines a new probability distribution. Likewise, if we 
superpose the corresponding quantum states, we obtain a new quantum pure state. 
But if we average over the coresponding pure state density operators, pj = 
there results a mixed state. The latter will be discussed in section 6. 

Now it is important to bring out a special form of state superposition. This is the 
Schmidt decomposition, 

M=E*ilV*>i® llfe>i, (3-18) 
i 

which exists for any bipartite state (see e.g. [6]). It must be emphasised that both 
factor states in the above tensor products may themselves correspond to several degrees 
of freedom, but the result is only proved if there are only two of them. The product states 
form a particular orthonormal basis in which to describe the state, so that the real, 
non-negative coefficients, \j, satisfy J2j <\j — 1- The state is entangled, unless \j = Sij. 
The Schmidt decomposition is often employed for the description of entangled states in 
finite Hilbert spaces. In this case, the number of nonzero eigenvlues, Xj, is a relevant 
quantifier of entanglement, known as the Schmidt number. For infinite dimensional 
Hilbert spaces, there may be an infinite number of nonzero Schmidt coefficients. 



4. Semiclassical quantum states 



Consider a momentum eigenstate \p') for L = 1. In the momentum representation, this 
is just 

(P\p')=t(p'-P), (4-1) 

which is not in a good form for semiclassical extrapolation. For this purpose it is better 
to use the complementary representation, 

.Sp>(q)\ 



(q\p') 



exp 



iqp 



exp 



(4.2) 



The phase in this expression can be interpreted as the area between the classical curve 
(the straight line, p' = p) and the g-axis. There is also an arbitrary constant phase, 
which is established by the choice of the initial point for the integral, 

rq 



S(q) 



p(q) dq. 



(4.3) 



Consider now a general observable, K(p,q). Its eigenstates correspond classically 
to curves, 7, in phase space: K(p, q) = k. These may be viewed locally as (possibly 
multivalued) functions, Pj{q)- Then the simplest semiclassical approximation is 



(#) = Y. A M) ex P 



+ V; 



(4.4) 
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Figure 2. The caustic of the semiclassical approximation to (q\k) lies in the 
neighbourhood of the point q c , where the tangent to the classical curve is vertical. The 
projection of this region onto the p axis is nonsingular, leading to a good semiclassical 
approximation for (p\k). 



see e.g. [19J. The phases, Sj(q) are again obtained from (14. 3p . The extra constant 
phases, Vj, are known as Maslov indices [HI [19], but they will not be discussed here. 
The amplitudes, Aj(q), are defined purely in terms of the classical structure. They are 
finite wherever the vertical line, q = constant, intersects the classical curve transversely. 
Where this vertical line is tangent to the classical curve, such as q c in FigfSJ the amplitude 
diverges. These points where the semiclassical approximation breaks down are known as 
caustics. The different branches of the function Pj(q) are connected at caustic points. 

In the case of bound eigenstates of K, the curves, 7, are closed. Then the eigenvalues 
are approximately obtained by the Bohr-Sommerfeld quantization condition, 



The quality of the semiclassical approximation for both the states themselves and their 
eigenvalues improves for large quantum numbers, n. Ground states, including that of 
the harmonic oscillator, are badly described by these approximations. 

Even for large n, a closed curve, 7, must inevitably have at least a pair of caustics. 
The way around this is to switch to the p-representation. Then the vertical tangent at 
the caustic position, q c , shown in FigfJl would correspond to the state, (p\q c ), which is in 
a nice semiclassical form. This means that the local branch of the multivalued function, 
q(p), gives rise to a semiclassical approximation which is a superposition of terms of the 
form 




(4.5) 



{p\k) = B(p) exp i 



S(p) 



+ v 



(4.6) 



h 
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This allows us to define the correct semiclassical approximation in the g-representation 
through the caustic region by the Fourier transform 

<«!*>= (SSJS* /<MP«exp(f), (4.7) 

which leads to a more refined approximation in terms of Airy functions instead of 
exponentials. This is usually refered to as the Maslov method of dealing with caustics 
[2"U] (also discussed in [T9"]). 

Let us now consider a product state for L > 1. Then, 

/i a fmp'A (mp'A (n-p'\ (AQ , 

(q\p) =exp I — — 1 ...exp I — — 1 = exp I — — 1 (4.8) 

and we can generalize the definition of action, 

S(q)= f q p(q)-dq. (4.9) 

Jqo 

This does not depend on the choice of path between q and q, because p'(q) is a constant 
in this simple case. Hence, this function defines a Lagrangian surface, i.e. a surface 
such that 

£p-dq = 0, (4.10) 

for any (reducible) circuit [13] . 

In general, the product state will involve arbitrary eigenstates of L observables, 
K = K1K2...KL, each in its own Hilbert space: 

(q\k) = (q 1 \k 1 )...(q L \k L ). (4.11) 

The wave function will be a superposition of terms with the form 

(#) = II4(«) «P l -[S 1 (q 1 ) + ... + S L (q L )\, (4.12) 
1 n 

one term for each branch of the funtions, Pi{qi). 

Defining again S(q) as the above phase, it is seen to be independent of the 
order in which we progress along each segment (qoi,qi), while keeping the other 
integration variables constant: The definition (14.31) . now reinterpreted as a path integral, 
is independent of the path on the surface. Therefore, this more general surface, 
K(p> <?) = k, is also Lagrangian. 

If the surface is the product of L quantized circles (closed curves), it will be an L- 
torus, r. Each of the L irreducible circuits, 7/, must then satisfy the Bohr-Sommerfeld 
conditions, 

/ Pi ■ d qi = (n, + hh, (4.13) 

or some suitable generalization (see e. g. [19]). Notice that the line-integral here used 
is not restricted to plane sections of r, because all topologically equivalent circuits on a 
Lagrangian surface must have the same action. 
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Figure 3. Each point within the rectangular caustic of a two-dimensional product 
torus is the immage of four phase space points under p-projection. 

Let us now evolve the product state semi classically. The basic result, due to van 
Vleck [21] , can be reinterpreted as the statement that classical and quantum evolutions 
commute. In other words, we can evolve classically each curve, 7/, if there are no cross 
terms in the Hamiltonian, so that the different degrees of freedom are decoupled. Each 
evolved observable then corresponds to K t (xi, t) = Ki(xi(xi , t), 0) and we approximately 
reconstruct the classically evolved state from the evolved torus, T/(t), which is the 
product of the ■jiit) : K[(xi,t) = k\. 

Notice that this classical evolution of products of curves fits in to the general view 
concerning the evolution of product probability distributions in the previous section, 
by merely choosing fi(t;xi) = 5(Ki(xi,t) = k[) and running time backwards. The 
important distinction between classical and semiclassical evolution is that the latter 
contains interferences between the different branches of the evolving classical curve. 
Each representation exhibits these interferences in a different way. 

Just as cross terms containing products of the different variables in the Hamiltonian 
destroy the product form of a classical probability distribution, the classically evolved 
L-D surface corresponding to an original product state also ceases to be a product. 
However, the smoothness of the evolution implies that the topology of the surface must 
be preserved(be it plane, torus, or, in between: cylindrical) . Furthermore, the classical 
evolution, canonical transformation, and hence all reducible circuits on 

the evolved surface have zero action, i.e. r t still has the Lagrangian property, which 
allows us to define the path- independent action S(q), and the irreducible circuits of r t 
still satisfy the same Bohr- Sommerf eld conditions, to first order in fi. 
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Figure 4. Catastrophe theory establishes the generic form for the unfolding of the 
double fold caustic at each corner of the projection of a product torus as it evolves. 

Let us investigate further the case of two degrees of freedom. The separable torus, 
r = 71 (g> 72, can be pictured through the separate 71 and 72 curves. These coincide 
with sections of the 2-D torus by alternative 3-D planes. (The normal case for Poincare 
sections, see e.g. [19]) The 72 curve does not depend on the choice of the q\ = constant 
section. The separable torus projects as a rectangle onto position space (gi, q 2 ), as shown 
in FigEl Within this rectangle, there are four different branches of the torus, which 
project onto each position, q, corresponding to the combinations of the two branches of 
each circle. The caustics at the side of the rectangle are double fold lines. 

After a general canonical evolution, the sections of r are no longer equal for different 
choices of qi = constant (or q 2 = constant), though all the sections have the same area, 
Si (or 5*2). In some cases (to do with time-invariance of the Hamiltonian) the projection 
onto the g-plane will merely distort the rectangle, which will still have finite-angled 
corners connecting double fold lines. But in general, these corners, hyperbolic umbilic 
points, will unfold in the generic form specified by catastrophe theory, as shown in FigJH 
There are four possibilities for the topology of the unfolding of the rectangle, shown in 
Figj5j For L > 2, the projection of the L-torus onto the L-D g-plane will be a solid 
hypercube that will be distorted, or unfolded by the motion generated by a coupling 
Hamiltonian. (These geometries are reviewed in [19], but are more thoroughly discussed 
in [22].) 

The representations of quantum states in terms of orthogonal position, or, 
alternatively, momentum eigenstates are the best that we can do, because of Heisenberg's 
uncertainty principle. Semiclassically, this corresponds to viewing a Lagrangian surface 
through a set of Lagrangian planes that foliate phase space. We switch from the q- 
representation to the p-representation by means of a Fourier transform of (q\ip)- This 
corresponds classically to taking the Legendre transform of S(q) [13J. For L > 1, 
we may take the Fourier transform for a subset of the degrees of freedom. This 
corrresponds to using a classical description in terms of the alternative Lagrangian 
planes (p u ...,p h q l+1 , q t ). 
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(a) (b) (e) (d) 



Figure 5. The full topology of the full fold lines is not determined by catastrophe 
theory: Each of the above forms corresponds to a different symplectic evolution from 
an initial product torus. 

One way to achieve a full phase space description is to use the basis of coherent 
states [61], [2H [251 S [27], labeled by the phase space vector, 77 = (r] p , rj q ), 

m = © 1/4ex p ( - Yn {q ~ + lV i {q - (4 - 14) 

Even though the coherent state basis is overcomplete, the exact decomposition, 

|V) = - fd V \ V )( V ^), (4.15) 

71 J 

is unique. The coherent states are phase space translations of the ground state of the 
harmonic oscillator (with unit mass): 

(,|0> = Q^exp(-^). (4,6) 

These result from the action of the translation operator: 




(4.17) 



using the skew product (13.101) . If either r\ v = 0, or r] q = 0, we obtain the usual 
translation operators for momenta, or positions, respectively. The arbitrary phase due 
to noncommutation of p and q is here chosen in the most symmetric way, using the 
Baker-Hausdorf relation [26] . 

In quantum optics it is customary to switch to the basis of creation and anihilation 
operators (q ±ip)/y2h. In this context, the translation operator (14.171) depends on 
the complex chords (£ p ± i£ q )/\r2H and is called the displacement operator j6T] . The 
semiclassical limit for a complex phase space is not as transparent as the real theory 
treated here. However, it is quite feasible to effect phase space translations in an 
experimental optical context. [29|. 

The coherent state representation is not orthogonal and is overcomplete. The 
alternative, to be explored in the next section, is to work directly with operators: We 
represent operators in orthogonal operator bases in analogy to the way that quantum 
states are commonly decomposed. This allows us to work directly with the translation 
operators, without having to apply them to the ground state of the harmonic oscillator. 
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5. Operator representations and double phase space 

The linear operators, A, that act on the quantum Hilbert space form a vector space of 
their own: \A)). Defining the Hilbert-Schmidt product, 



should note the similarity between this dyadic basis, in the case of L = 1 with 

the basis of product states, \qi) <8> \q.2l- The substitution of a bra by a ket in the former, 
will in most cases imply no more than complex conjugation. 

Thus we may relate the vector space of quantum operators to a double Hilbert space 
with respect to that of quantum states. Since we have explored the correspondence of 
the state-Hilbert space with classical phase space, it is now natural to relate the double 
Hilbert space to a double phase space : X = x_ x x + (se e.g. [31]). The operator \Q)) 
should then correspond to the Lagrangian plane, Q = constant in the double phase 
space. This does hold, within a minor adaptation, analogous to the use of the adjoint 
operator in the definition of the Hilbert-Schmidt product. That is, we should define 
Q = 9+)) but P = (—p_,p + ) as coordinates of the double phase space X = (P, Q). 

A good reason for this is that then we include among the set of Lagrangian 
surfaces in double phase space all the canonical transformations in single phase space, 
x_ — > x + = C(x_). This also transports closed curves, 7_ — > j + , so that we may rewrite 
the definition of a canonical transformation ( 13.61) as 



where F = (7_,7 + ). Thus we may consider j± as projections of the curve F defined on 
the (2L)-dimensional surface, Ac, which specifies the canonical transformation, within 
the (4L)-dimensional double phase space, X = (P,Q). 

It is worthwhile to consider the richness of structures in double phase space. On 
the one hand, a canonical transformation defines a Lagrangian surface as x_)-(x_), a 
one-to-one function. On the other hand, the product of a Lagrangian surface, A_ in x_ 
with another surface A + in x + , A = A_ ® A + , is also Lagrangian in double phase space, 
but projects singularly onto either of the factor spaces. In the case that both sufaces 
are tori, we obtain a double phase space torus, r = r_ ® r + , as if we had doubled the 
number of degrees of freedom. (All Lagrangian surfaces will hereon be labled r, even 
when they are not necessarily a torus; in the case that L — 1, r is just a closed curve, 
7.) If L = 1, it will be a 2-D product torus, with the only difference that p- — > —p- 
in the present construction. If each Lagrangian surface corresponds to a state, i.e. 
and 1^+)) then we represent \^/)) = \ip + ){ip_\ in the \Q)) representation as 






(5.3) 




(5.4) 
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Therefore, the semiclassical approximation is just a superposition of terms of the form 
((Q\n = Aj(Q)exp[iSj(Q)/H], (5.5) 



with 



and 



A J (Q) = A j ^y A j+ (q + ) (5.6) 



Sj(Q) = £ Pj(Q')-dQ'. (5.7) 

Again this is in strict analogy to the construction of semiclassical product states of 
higher degrees of freedom. Note that the projection of the double Lagrangian torus 
onto P, or Q is just the rectangle discussed previously for product states, whereas the 
projections onto the planes, x_ and x + are specially singular. 

The semiclassical approximation for a unitary operator, U, that corresponds to a 
canonical transformation, C : X- — > x+, has exactly the same form, i.e a superposition 

((Q\U)) = (q+\U\q-) = Uj{Q) e W [iSj(Q)/K], (5.8) 

for each branch of the function Pj(Q) defined by the Lagrangian surface in double phase 
space. Note that the situation with respect to projection singularities is now reversed, 
as compared to \^))- The fact that the projections of the Lagrangian surface, A c , onto 
either x_, or x + are both nonsingular in no way guarantees that the projections onto 
the P, or the Q Lagrangian planes will be likewise free of caustics. 

Conversely, any function, S(Q), is, at least locally, the generating function of a 
canonical transformation through the implicit equations: 



Here we recognize the standard generating functions S(q-,q + ) in e.g. Goldstein [16J. 
If S(Q) is quadratic, then these implicit equations will be linear, so that the explicit 
transformation will result from a matrix inversion (if it is nonsingular). There will be 
a single branch in S(Q) for such a symplectic transformation and it turns out that the 
semiclassical approximation is exact in this case. 

The well known alternatives to these generating functions are usually obtained by 
Legendre transforms. However, we can consider the 7r/2 rotation, q + — > p + ,p + — > —q+, 
times the identity in x_, as an example of canonical transformation in double phase 
space: X — > X'. Then Q' = (q_,p + ) is also a good Lagrangian plane which can be 
used as the new coordinate plane for the description of Ac- In the new coordinates the 
implicit equations for the canonical transformation are just 

|| = P W ), or ^ = £ = _ p _. (5,0) 

The correspondance with a semiclassical state, 

«Q'|tf » = A'j{Q') expiiS'jiQ^/H], (5.11) 
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will be exact in the case of a symplectic transformation. Note that \Q')) is a first example 
of an operator basis that corresponds to a set of parallel Lagrangian planes in double 
phase space, which, nonetheless, have internal coordinates that can be identified with a 
phase space on its own. 

The crucial step is now to explore other kinds of canonical transformations in double 
phase space [32]. In particular, 



Here, the J symplectic matrix in single phase space is essential to canonize what would 
be just a 7r/4 rotation. It accounts for the change of sign in the p_ coordinate. We will 
here have to bare the disconfortable situation that the canonical coordinate in double 
phase space is y, but the geometrically meaningful variable in single phase space is £, 
the trajectory chord, in the case of continuous evolution. The coordinate x will be 
referred to as the centre. 

If we consider the horizontal Lagrangian planes y = constant, each is identified 
with a uniform classical translation. Thus, we have departed from coordinate planes 
corresponding to dyadic operators to those planes in double phase that describe 
canonical transformations and hence correspond to unitary transformations. In this 
case, X— = x_ ~\~ £ are the group of phase space translations, which include the 

the identity, i. e. the identity plane is defined as £ = 0. 

On the other hand, the vertical plane, x = 0, defines the canonical reflection 
through the origin, x_ — > x + = —X- (or inversion), since all the chords for this 
transformation are centred on the origin. Other vertical planes specify reflections 
through other points, X- — > x + = — (x_ — 2x). The reflections do not form a group on 
their own (no identity), but together with the translations they form the afRne group 
of geometry [33J . 

Since there is an exact correspondence between linear canonical transformations 
and unitary transformations, each plane y = constant corresponds precisely to the 
translation operator, T^, previously defined as ( 14.171) . Notice that this was written with 
a phase that is a skew product involving £, but we could also use Tg = exp[iy ■ x/h]. In 
terms of the previous dyadic \Q)) basis, this is expressed as 



a symmetrized Fourier transform (see e. g. [32]) • 

Just as a 7r/2 rotation in single phase space, q — > p, p — > — q, corresponds to a Fourier 
transform, so the transformation between horizontal and vertical planes in double phase 
space is also achieved by a full Fourier transform (except for an annoying factor of 2 L ): 




(5.12) 




(5.13) 




(5.14) 



In terms of the dyadic \Q)) basis, we have 




(5.15) 
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A = I TK^Z m % (5.19) 



the complementary symmetrized Fourier transform to (15. 13j) . 

We are now free to switch from the usual (position) dyadic basis to the unitary 
operator basis, |y)) = T^. 

{(y\A)) = trf -tA = A(t), (5.16) 

where we use T_£ = T^. A{£) is the chord representation of the operator A. (Also 
referred to as the chord symbol). To verify that this is indeed the expansion coefficient 
for an arbitrary operator in the basis of translation chords, we use 

trf € =(27r«)^(0 = «y|/» (5.17) 

(note the double phase space analogy with (p'\(p = 0)) = 5(p') ), as well as the quantum 
version of the group of translations: 

44=4+6exp[^kiA6] (5.18) 
(see e.g. [S]). Then, the expansion 

_^ 

(2vrfi) 

leads to 

= tr / A «' )f = / m exp [ k = m - (5 - 20) 

The chord representation is thus a second example of a representation of operators in 
terms of an operator basis that can be identified uniquely to a phase space. Indeed, each 
chord corresponds to a Lagrangian surface in double phase space and hence a particular 
uniform translation in single phase space. 

The next representation will be based on phase space reflections, i? x . But first it is 
worthwhile to examine some characteristcs of these operators. Unlike the translations, 
they do not form a group on their own, though they combine with the latter to form 
the affine group. The products are [31] 

i? x T ? = exp[-^x A i] R\^ /2 , (5.21) 
n 

f^R x = exp[~x A f ] R xH/2 (5.22) 

and 

_ _ 2i 

-Rx 2 -R Xl = exp[— xi A x 2 ] T 2(X2 _ X1 ). (5.23) 
n 

Except for the phases, these are just the classical relations. The last one is specially 
interesting. Note that i? x = /, the identity, hence the (degenerate) eigenvalues of _R X 
must be either +1, or —1. Therefore these operators are Hermitian, as well as unitary. 

Are they true observables? Consider the effect of Rq on the eigenstates of the 
harmonic oscillator. Taking q — > —q and p — > —p, leads to a change of sign for all 
the odd states, while preserving the even states. In other words, the latter are just the 
(+l)-eigenstates, while the odd states are (— l)-eigenstates. Though it is hard to imagine 
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measuring the parity of a particle, we saw in section 2 that the parity decomposition of 
even a classical wave can certainly be effected. Measurements of the eigenvalues of this 
non-mechanical observable are currently performed for single photons in optical cavities 
|30j . It is true that these measurements are performed on a mode of the electromagnetic 
field raher than a particle, but it only makes sense to discuss the parity within a specific 
mode if it is quantized. 

Reflection operators are very strange observables as far as phase space 
correspondence is concerned. It was discussed in section 4 that usual observables 
correspond to smooth phase space functions and their eigenvalues correspond to 
level curves if L = 1. This is just not the case of reflection operators with their 
infinitely degenerate, ±1, eigenvalues. In their dual role as both unitary and Hermitian 
(observable) operators, reflections are almost schizophrenic: They are perfectly ordinary 
unitary operators, corresponding to Lagrangian planes in double phase space, but they 
do not correspond to any smooth classical function in phase space, as expected of a 
mechanical observable. 

This should furnish sufficient motivation to investigate the representation of 
arbitrary operators in terms of reflection centres. The assumption that 

A= f dxA(x) 2 L i? x , (5.24) 



leads to 

«x|A» = tr (2 L R X )A = trj ^^A(x')(2 L R x )(2 L R x ,) = A(x). (5.25) 

This is the Weyl representation of the operator A (also known as the Weyl symbol). 
Once again we use half the coordinates of double phase space, within a Lagrangian plane 
that is a phase space on its own, to describe a quantum operator. This perception that 
we are really dealing with different phase spaces for each operator representation was 
clearly stated in the excellent review by Balazs and Jennings (36] . What was lacking was 
merely the identification of each of these different phase spaces with a specific foliation 
of Lagrangian planes in double phase space. 

As far as unitary operators, U, are concerned, the semiclassical limit of the 
representations, either in terms of centres, or chords has exactly the same form as 
for any other Lagrangian basis. For instance, the Weyl symbol will be a superposition 
of terms, such as 

U(x) = A(x) exp[iS(x)/fi], (5.26) 
in terms of the centre action, defined as 

5(x) = / X y(x') • dx' = r^(x') A dx'. (5.27) 
Jo Jo 

For symplectic transformations, the Lagrangian surface is a plane and so there is only 
a single branch of the action function 5*(x), which is quadratic. Then (15. 26ft is an exact 
representation of the corresponding quantum metaplectic transformation. However, 
in the general nonlinear case there may be caustics in the projection of the Lagrangian 
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y (x') surface onto the x-plane. Recall that this is just the plane that defines the identity 
operator, / (corresponding to y = 0, or S = 0). 

For the canonical transformation generated by a Hamiltonian, H(x), it turns out 
that the generating function has the limit [M] 

5 e (x, t = e) -> -eff(x) + C(e 3 ). (5.28) 

There are no caustics for small times in the centre representation, since the 
corresponding Lagrangian surface is nearly horizontal. 

The smooth real Hamiltonian itself can be equated to the Weyl symbol for the 
corresponding operator, H, within semiclassically small ordering terms. This is the case 
of the Weyl representation for any observable that corresponds classically to a smooth 
classical function of the points in phase space [TTj . Since we can always consider classical 
observables as infinitesimal generators of motion through Hamilton's equations, it is 
appropriate to picture them as functions on the y = plane, so that the Hamiltonian 
vectors form a field on this plane that indicates which way it will evolve. In contrast, 
the chord symbol for these smooth mechanical observables is not at all smooth. This 
is because the chord and centre symbols are related to each other through the Fourier 
transform, 

A(0 = / dx exp (-l£ A x)A(x), (5.29) 

just as the translation and reflection operators themselves in (15. 14ft . This Fourier 
transform takes the symbol for the identity, /(x) = 1, into = and a Taylor 
series in x into a series of derivatives of (^-functions. However, we shall see in the next 
section that the chord representation of density operators have very useful properties. 

It is fitting to consider here another feature which distinguishes the reflection 
operators from mechanical observables. Far from being represented by a smooth phase 
space function, their centre representation is just 

i? x (x) = 2- i 5(x / -x). (5.30) 

These singular functions cannot be interpreted as corresponding to classical states (i.e. 
individual phase space points) because the _R X have the eigenvalue —1, so they are not 
density operators. 

Probably the first to remark on the general structure of translations and reflections 
underlying the Weyl and the chord representations were Grossmann and Huguenin [35J. 
There exists an exact correspondence, between these operators of the affine quantum 
group, together with the unitary operators of the metaplectic group with the classical 
transformations of the inhomogeneous symplectic group [T7]. In other words, all linear 
canonical transformations, including reflections and translations, are exactly matched by 
quantum unitary transformations. Thus, the unitary transformation, Uc, corresponding 
to x — > x' — Cx, where C is a symplectic matrix, takes 

R X ^R' X = U^RjJc = R*' and T s -> f£ = W c f^U c = (5.31) 
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This has the consequence that both the centre and the chord representations are 
invariant with respect to metaplectic transformations, because the transformed operator, 
A — > A' is represented by 

A'(x) = tr U C A U ] C R^ = tr A U ] c Rjj c = tr A R' x = A(x') (5.32) 

and, likewise, A'[£) = A(£'). 

This section is concluded with some general formulae concerning these 
representations. For the trace of an operator, we have the alternative forms: 

tr A = tr I A = ((T 5=0 |A» = A(£ = 0) = J dxA(x). (5.33) 

The adjoint operator, A\ is represented by 

At(x) = [A(x)]*, or At(0 = [A(-0]*, (5.34) 

where * denotes complex conjugation. Thus, if A is Hermitian, A(x) is real, though 
may well be complex. The Weyl or chord symbols for products of operators is not 
at all obvious (see e.g. |34j), but 

tr AA = / (^ A ^-v = / p^w* <*>■ < 5 - 35 > 

6. The Wigner function and the chord function 

It is customary to alter the normalization of the centre and the chord symbols for the 
density operator, p, so as to define 

, tt and y It) = -: 

respectively the Wigner function and the chord function. Combining with the 
general definition of the Weyl representation (I5.24[) and the expression (I5.15P for the 
reflection operator, we obtain the original definition of W(x), proposed by Wigner 
[37]. In both cases of (jfi.ip . the representation of the trace of a product leads to the 
expectation of any observable, A, as 

(A) = J dx W(x) A(x) = J d£ X (-0 (6.2) 

The first integral is more interesting because A(x) is at least semiclassically close to 
the classical variable, which tempts us to identify the Wigner function with a nearly 
classical probability distribution. However, we will see bellow that VT(x), though real 
and normalized so 

dxW(x)=l, (6.3) 

may well take on negative values. 

The chord function behaves like a classical characteristic function, in as much as 
the moments are 

(,») = trcf p = (ihr ^-(27rh) L X (0 t _ n (6.4) 



= 7^1 -d X (0 = (6.1) 
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and 

(p n > =trp n p= (-i«) a |lr( 2 ^) L X(0| £=0 • (6-5) 
Taking the zero'th moment, we obtain the normalization, 

1 = (2nh) L X (0), (6.6) 

because tr p = p(£ = 0) = 1. 

Shifting the phase space origin to (x) = ((p), (?)), we can define the Schrodinger 
covariance matrix [38] just as its classical counterpart ( 13. 7\\ . with <5p 2 = (jS 2 ), 
#g 2 = (g 2 ) and {Spq) 2 = ({pq + qp)/2). It is then obvious that the expansion of the 
chord function at the origin is just 

X (0 = -£K£ + ... (6.7) 

and we can interpret the uncertainty, 



A K = Vdet K, (6.8) 

as proportional to the volume of the ellipsoid: £ K £ = 1. Evidently, this volume is 

invariant with respect to symplectic transformations, so that Ak is a symplectically 
invariant measure of the uncertainty of the state. 
The projection of the Wigner function, 

dp W(p,q) = Pr(g) (6.9) 



is a true probability for position measurements [37]. Furthermore, the invariance of the 
chord and the centre representations with respect to symplectic transformations then 
guarantees that the projection of the Wigner function along any set of Lagrangian planes 
p' supplies the probability distribution for the conjugate variable q'. In particular, the 
probability Pr(p) results from the projection of H^(p, q) with respect to q. All these 
planes are Lagrangian, so it follows that the projection of the Wigner function onto 
any Lagrangian plane in phase space is a probability distribution for the corresponding 
variable. It may appear somewhat contrived, as far as measurement is concerned, to 
consider general linear combinations of position and momentum. However, it should be 
recalled that these observables will evolve from an initial position for the motion driven 
by any quadratic Hamiltonian, even including free motion through a laboratory. The 
reconstruction of the Wigner function from a suitable set of these marginal distributions 
is known as quantum tomography. This is achieved through the Radon transform (see 
e. g. [39]). 

It is equally remarkable, but less well known, that the characteristic function 
corresponding to the marginal probability distribution for positions is obtained by 
merely taking a section of the chord function: 

J dq Pr(q) exp {--^ ■ q) = {2nh) L X (0,^). (6.10) 

Since, the chord function is also symplectically invariant, it follows that the characteristic 
functions for all the probability distributions, which result from Wigner projections onto 
Lagrangian planes, are equal to the corresponding sections of the chord function. 



Entanglement in phase space 



27 



So far we have emphasised the seemingly classical aspects of the Wigner function. 
However, it must be remembered that the Weyl representation is defined in terms of a 
very anomolous observable, as far as classical correspondence is concerned. In order to 
reveal the full quantum nature of the Wigner function, let us divide the Hilbert space 
of quantum states into even and odd subspaces for a given reflection operator, P x . This 
is achieved through the projection operator introduced by Grossmann [ID] and Royer 

P£ = ±(l±R x ), (6.11) 

so that, in its turn, we can express each reflection operator as the superposition of this 
pair of projections onto the even and the odd subspaces: 

P X = P*-P X . (g 12 ) 

But 

trpP£ = Pr£ (6.13) 
is just the probability of measuring P x to have the eigenvalue ±1, so it follows that [JT] 



We thus find that the Wigner function does not admit the interpretation as a probability 
distribution in phase space, because it can certainly be negative. Even so, it is a simple 
linear function of a distribution of probabilities of positive eigenvalues for all possible 
reflection measurements. Its maximum possible value (ttH)~ l is attained for any point, 
x, such that P* p = p, whereas the commutation of the density operator with P x 
specifies a phase space point where W(x) = — {ixti)~ L . 

Let us now investigate the effect of reflections and translations on a density operator. 
Evidently, the centre and chord representations are specially suitable for this purpose. 
In the case of a phase space translation by the vector, 77, i.e. p v = T v pT- v , the respective 
Wigner and chord functions become 

W„(x) = W(x - 77) and x v (0 = ( 6 - 15 ) 

which shows that, unlike the Wigner function, the chord function is not generally real. 
The sensitivity of a state to translations is described by the phase space correlations 
of a given density operator, defined as 



C(f) =trp%pf\ = (27ih) L J dxW(x) W(x-0 = {27Th) L J dn e iv ^ /h \ X {v)\ 2 ■ (6.16) 

^From the reciprocal relation that supplies the intensity of the chord function as the 
Fourier transform of these correlations and the normalization condition (16. 6ft . we see 
that 

*d V C(rj) = (2nh) 3L \ X (v = 0)| 2 = (2vrfi) L . (6.17) 

So, even though these correlations are defined in terms of classical translatons in phase 
space, they are purely quantum and disapear in the classical limit. However, if we 
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fix h and adopt this constant as our phase space scale, then we can picture C(£) as a 
classical-like phase space distribution for which the characteristic function is just |x(OI 2 - 
Specializing to the case of a pure state, p = \ip)(ip\, we find that 

(V>|TdV> = (2vrfi) L X (-0 , (6.18) 
so that the phase space correlations take the form [42] 

C(0 = \(mm 2 = (2nH) 2L \ X (0\ 2 - (6.19) 

Thus, for instance, in the case that £ = (0,£ 9 ), the phase space correlations, C(£), 
are just the usual spacial correlations inferred from neutron scattering experiments. 
Nonetheless, we must be carefull to distinguish between phase space correlations and the 
correlations between the quantum measurements of observables defined on the different 
components of a bipartite system, such as enter the CHSH inequality. For a pure state, 
(I6.19P is the square modulus of the expectation for a translation, which is not a quantum 
observable. However, (16.161) defines the phase space correlation in the same way as for 
a classical distribution. 

The chord function always assumes its maximum value l/(27rfi) L at the origin. But 
also an average of overlaps cannot exceed one, so x(0) is the maximum even for mixed 
states. As for the correlations, we always have 

tip 2 = (2nH) L J dx [H/(x)] 2 = (2nH) L J d£ | X (0| 2 = C{0), (6.20) 

being that tr p 2 = 1 for pure states. But consider a mixture of orthogonal states, 

p= ^Pr(n) \n)(n\, (6.21) 

3 

then the purity 

C(0) = ^Pr(n) 2 < 1. (6.22) 

3 

Another form in which this quantity appears is the linear entropy : 1 — tr p 2 . This 
may be considered as a first order expansion of the von Neumann entropy: 

- tr p In p = - J dx W(x) ln[(27rfi) L W(x)] , (6.23) 

a quantum version of the classical Shannon entropy. In [42J the correlations were 
normalized by the purity so as to be always unity a the origin, but it is convenient 
to include this quantity as a special case of the correlations. 

General invariance with respect to Fourier transformation characterizes the 
correlation in the case of pure states. Indeed, inserting the above expression in the 
definition of the phase space correlation (16.161) . we obtain [12] 

This is a remarkable property of all pure states and is in no way restricted by special 
symmetry properties that will be shown to relate certain Wigner functions to their 
respective chord functions. An immediate consequence is that oscillations of the phase 



Entanglement in phase space 



29 



space correlation of a pure state involving a large displacement, £, are necessarily bound 
to small ripples on the scale, in the direction, J£. Of course, these small scale 

oscillations of the phase space correlations, which have been attractively described as 
subplanckian [13], show up in the pure state Wigner function because of (16.161) . 

The Fourier invariance condition (16. 24ft includes as a special case the more familiar 
one obtained by tracing over the full pure state condition p 2 = p. It follows that the 
difference of both sides of (16.241) for each chord £ generalizes (16.221) as a measure of the 
degree of purity of a state. All the same, the loss of the phase information in C(£), 
but contained in the chord function, would seem to imply that these are necessary 
conditions, whereas the full sufficient condition of purity is p 2 = p, which is expressed 
in the chord representation as [34J 



d?7 X (v) X(£ - V) e i?Ar?/2fi = J d V x^(v) x(£ ~ v) = x(0 , (6.25) 

with x$/2(v) defined by (I6.15p . However, the particular condition C(0) = 1 is indeed 
a suficient condition, because, for any mixture of pure states, p = J2 n P T ( n )\4>n)(4>n\, we 
obtain 

tr f = 1 - J2 Pr(n)Pr(n')[l - \(ip n \M\ 2 )- (6.26) 

A single phase space point does not correspond to any pure state in Hilbert space. 
The only pure states that are classical-like, i.e. have positive Wigner functions, are 
either coherent states, or their immage by a symplectic transformation [4~4"t |4"5]. 

Let us now consider the effect of measuring a general phase space reflection, R x . 
The density operator, p, will be projected by P±, defined by (16.111) . onto either the even, 
or odd subspace for this particular relection: 

P x n P x 

pi = ■ (6 - 27) 

tr p 

The Weyl symbol for p?| coincides with the symmetric Wigner function, W^(x'), 
within a normalization factor, so that we obtain 



^(x<) = W-trfl- (^) PI = ^ ± . (6 . 28) 

It follows that the Wigner function and the chord function for a reflection symmetric 
density operator are trivially related. Shifting the origin of phase space to the symmetry 
point leads to [12] 

W^(x) = ± 2 L X ±(-2x) . (6.29) 

Thus, all Wigner functions for density operators that commute with a reflection 
symmetry attain the largest amplitude at the symmetry point, but this will be negative 
in the case of odd symmetry. 

% For distributions, Pr(n), over eigenstates \n) of an observable with discrete spectrum, the condition 
Pr(n) 2 = Pr(n) also singles out a pure state, Pr(n) = Sn,m-> but this condition is not generalizable to a 
continuous spectrum. 
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Let us consider some standard examples of Wigner and chord functions. All the 
following cases are related to eigenstates of a harmonic oscillator with one degree of 
freedom and unit mass. 

i) Coherent states: the Wigner function is just a Gaussian centered on 77, 



W„(x) 



whereas, 



7rfi 



x,(0 



2nh 



cxp 



exp 



H 



1 



yP ~ Vp) 



u=l 



irh 



exp 





2nh 



(6.30) 



. (6.31) 



So, any translation of the coherent state merely alters the phase of the Gausssian 
chord function that sits on the origin. The coherent states, or more generally all 
equivalent Gaussian states obtained from them by symplectic transformations, are the 
only examples of pure states for which the Wigner function is nowhere negative [44]. 
This is one of the reasons why these are sometimes considered to be the most classical 
of pure quantum states. Since the projection of a Gaussian is also a Gaussian, the 
measurement of position, or any other Lagrangian phase space coordinate, does not 
display interference fringes. The fact that the uncertainty, A = 5p5q = H, is minimal 
allows us to interpret them as quantum phase space points. 

(ii) A superposition of a pair of coherent states, \rj) ± | — rj) is sometimes known as 
a Schrddinger cat state. Its Wigner function is El 



W±(x) 



e -(x-,)Vfi + e - 



± 2e- 



2 /H 



cos -x A rj 

h 



(6.32) 



2 7 rfi(l±e" r ' 2 /fi) 

It consists of two classical gaussians centred on ±i] and an interference pattern with a 
gaussian envelope centred on their midpoint. The frequency of this oscillation increases 
with the separation \2rj\. In fig]6]the displacements, ±77 have been chosen as (±3, ±3). 
The phase of the pair of coherent states merely shifts the phase of the interference fringes, 
so that the midpoint is an absolute maximum for W + (x) and an absolute minimum for 
W_(x). It might be supposed that, for small rj — > 0, we would have VF+(x) > for 
all x, but it is easy to verify that there are very shallow negative regions far removed 
from the classical superposed Gaussians, in agreement with [44] 145] . The interference 
pattern of the Wigner function does not survive the projection orthogonal to rj: In this 
direction, the interference disappears to produce a purely classical pattern. Conversely, 
the projection along 77 is marked by interference fringes. 
For the chord function, 



x±(0 



1 



-{ii/2-nf/h 



+ e 



± 2e 



cos -£ A rj 
H 



(6.33) 



this same configuration has to be reinterpreted. Now the local phase space correlations 
of the individual coherent states, as in (i), are placed in the neighbourhood of the origin, 
where they interfere, while their cross-correlation generates new Gaussians centred on 



Here and below we set lo = 1. 
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Figure 6. The Wigner function for the Schrodinger cat state displays a pair of classical 
Gaussians, one for each coherent state, and a third Gaussian modulated by interference 
fringes halfway between them. The chord function is a mere rescaling of the Wigner 
function if the midpoint lies on the origin. 
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the separation vectors ±2i]. The general case of coherent states \rjx) and (772) merely leads 
to Gaussians centred on ±(771—772) with addition of the phase factor expfi^+T^) A£/2fi]. 

Recalling that the phase space correlations of a pure state are just the square 
modulus of the chord function, we can immediately verify the general relation betweeen 
large and small scale structures ( 16. 24ft in the case of Schrodinger cat states. Indeed, the 
spacial freequency of the oscilations of the chord function increases directly with the 
separation of the pair of coherent states. 

The particular superpositions of coherent states, |+) and |— ), are respectively even 
and odd eigenstates of the parity operator R , i.e. reflection about the origin. Therefore, 
they are the two possible states that could be produced by a parity measurement effected 
on the single coherent state \r]). Thus, the parity measurement would generate a sizable 
probability of finding a particle near x = —rj, even though this were most unlikely before 
the measurement. 

The states |±) are orthogonal, even though the coherent states, \rj) and | — rj), 
are not. It is true that such a pair of coherent states will be nearly orthogonal if 
rj is large enough and thus considered to form a qubit. Within this approximation, 
the symmetrical states would then be a mere unitary transformation of a single qubit. 
However, no approximation is needed in this process of carving a qubit from an infinite 
dimensional system, if we use |±) as the original basis states. We would then consider 
a common garden coherent state to be the superposition of a symmetrical pair of 
Schrodinger cats. (Is there some approximation involved?) Indeed, this generation 
of a qubit by a reflection is not limited to coherent states, but could in principle be 
realized for any unsymmetrical initial state. 

(iii) Fock states, \n), i.e., the excited states of the harmonic oscillator, also have 
reflection symmetry with respect to the origin. Thus, from the exact Wigner function, 
first derived by Gronewold [47J, 

W B (x) = ^"* 2 / fi L n ^pj , (6.34) 

where L n is a Laguerre polynomial, we obtain the chord function 



e 



-5 2 /4fi 



*■<« = -srMsJ ■ (6 ' 35) 

It is interesting to note that the symmetry centre, which produces the maximum 
amplitude of the Wigner function, is nowhere near the classical manifold with energy 
E n = (n + |) Hlj. However, this point lies in a region of narrow oscillations, so that it 
does not affect the average of smooth observables. FigJT] shows the Wigner function for 
the Fock state with n = 2; the origin is a maximum because of the posititive parity. 
The unfolding of this peak for nonsymmetric Wigner functions is discussed in section 
10. 

The Wigner function exhibits the interference fringes for the measurement of any 
variable ap + bq. In the case of the Fock state, these are always present. A simple way 
to see this is that any direction for the projection will be somewhere tangent to each of 




Figure 7. The Wigner function for the n = 2 Fock state. The classical Bohr-quantizcd 
circle lies just outside the maximum of the outer fringe. The chord function is just a 
rescaling of the Wigner function. 
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the continuous curves that form the Wigner function fringes. These regions dominate 
the projection. This example thus illustrates the necessity for the Wigner function to 
have negative regions: This is the only way that interference can result from a mere 
projection in phase space. 

The Fock states are an example of a complete parity basis, which is even or odd 
according to the state label, n. Hence, if a pure state, is specified in this basis, then 

W*{0) = \j2l\(^)\ 2 ~ l((2n + 1MI 2 ]. (6-36) 

Such a decomposition can in principle be achieved for general arguments of the Wigner 
function, but then it is necessary to translate the whole Fock state basis instead of just 
the ground state, as in the definition of coherent states. If we similarly translate the 
Hamiltonian, it will commute with i? x instead of commuting with R . The eigenstates 
of all such Hamiltonians will form a good odd-even basis. The difficulty with defining 
a semiclassical correspondence for both these classes of eigenstates is that the odd and 
even Bohr-quantized curves approach each other without a limit as H —>■ 0. 

All the above examples are singled out by some point of reflection symmetry, which 
needs to be chosen as the origin for the chord function to be real. The chord function 
must assume its maximum value l/(27rfi) L at the origin, whatever the symmetry, because 
of normalization. The Wigner amplitude, |W(x)|, need not have such a prominent peak 
in general. However we shall see in section 10 that the large scale features of the 
semiclassical forms of the Wigner function and the chord function maintain a mutual 
correspondence, even in the absence of a reflection symmetry. 

It is important to note that the commutation, [p, Rq] = 0, guarantees that W(x) 
is a symmetric function with respect to (classical) reflection at the origin. This is 
a consequence of the fact that, if pr = R p R, then W^x) = W(R(yi)), the classical 
reflection of the argument. However, it is the maximum (or minimum) value at the origin 
which guarantees that the density operator is pure with respect to parity, i.e. it is either 
p + , or p_. Indeed, even though a mixture of an even density and an odd density (i.e. 
p = c + p + + c_p_) will trivially satisfy [p, R ] = 0, we see that W(0) = (c+ — c_)/(7rfi) L 
will not be maximal. 

Perhaps the converse property is of even more interest: If a Wigner function is 
symmetric about the origin, but |VT(0)| < (ttH)~ l , then the state must be a mixture. 
After all, if it is a mixture of parities, it cannot be a pure state. It is only if the mixture 
is restricted to states of the same parity that it will not be detected by W(0). There 
are many measures of degrees of mixedness, or impurity, but it is specially nice to be 
able to spot this property by a mere glance at the Wigner function. Furthermore, the 
Wigner function, i.e. the parity decomposition, is a measurable property [281 1291 130]. 

The Wigner function may be considered as a field of probabilities for parity 
decompositions in phase space. Each reflection separates the infinite dimensional Hilbert 
space into a pair of orthogonal components. If we just consider a single reflection, this 
goes a long way to reducing the Hilbert space to that of a single qubit, a two state system. 
No matter how classical the appearance of the Wigner function (i.e. it may be smooth 
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and positive) it is always fully quantum as far as parity measurements are concerned. 
The situation is quite different, for instance, for position measurements. Then there is 
an important difference between the Wigner function for a pure Schrodinger cat state 
and a mixture of cats with different phases. This is revealed by the fine interference 
fringes between the two classical regions, but is even more clearly displayed by the pair 
of correlation peaks far from the origin of the chord function. The relation between 
these features and entanglement is discussed in section 8. 

There is a vast litterature concerning the Wigner function. Only a few topics have 
been mentioned here and it has been necessary to leave out even as relevant a topic as 
quantum tomography. The adaptation of the Wigner function for finite Hilbert spaces 
is of special relevance for quantum computing and quantum information theory. Then 
the rule that each quantum state corresponds to a volume of (27rfi) L in classical phase 
space, restricts the overall phase space volume. Thus, one must first face a choice of the 
topology in which to compactify phase space. It turns out that the simplest choice is a 
torus, though single qubits are more naturally displayed on a Block sphere. In spite of the 
intrinsic interest in many of the aspects of finite space Wigner functions [i8j H5]; there 
is no overall agreement on the choice of Wigner function propertise to emphasise. Not 
all formalisms lead to a corresponding natural definition of a conjugate chord function 
as in [50], nor is there an overall preocupation with invariance with respect to those 
symplectic transformations which preserve the torus topology of phase space [51J. A 
final difficulty concerns the appearance of ghost immages and dimensionally dependent 
features [52] . 

So far nothing has been said of an alternative phase space representation, the Husimi 
fuction [53l 154]. Defined in terms of coherent states \rj) as 



because of (15.351) . The lack of purity of a state can be described in terms of the Wehrl 
entropy, 



According to Wehrl's inequality [56] (see also [57]), the Wehrl entropy is always bounded 
from below by the von Neumann entropy (16.231) . For more recent developments 
concerning the Wehrl entropy, see e.g. [58] . 

The Husimi function is most appropriate for the study of quantum chaos, because 
it highlights the classical region. But such a downplay of the quantum interferences, 
achieved by coarse graining the Wigner function, is not what one would ordinarily 
seek in quantum information theory. In a way, this is just the opposite of the chord 
function, which squashes all classical structure to the neighbourhood of the origin, so as 
to display the purely quantum coherences. It is remarkable that both these antithetical 



Ph{v) = (v\p\v) = tep\r])(v\, 
it can be interpreted as a smoothed Wigner function, 



(6.37) 




(6.38) 




(6.39) 
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representations are intimately related to the translation operators, since the Husimi 
ruction for a pure state, p = can be rewritten as 

Ph(v) = Mv)\ 2 = MT v \(v = 0))\ 2 - (6-40) 

Hence, the basic difference with respect to C v in (16.191) is the exchange of |0), the 
Gaussian ground state of the harmonic oscillator, for \%p) itself. 

A further comment is that the quantum interferences are displayed by the isolated 
zeroes of Husimi functions [55J, in the case that L = 1. A uniform distribution of zeroes 
has been used to characterize the eigenstates of classically chaotic systems. Even though 
this is of great theoretical interest, these zeroes are ususally located in regions where 
the Husimi function is already tiny, so that they may be very hard to compute. For 
instance, in the case of the cat state (16.321) with small n, they are found in the shallow 
negative regions where |W + (x)| is exponentially small. 



7. The partial trace: sections and projections 

Recall that the representation of operators, A = \A)), in a given basis, such as ((Q\A)), 
corresponds to the folliation of the double phase space, X = (P, Q) , by a set of 
Lagrangian planes, Q = constant. Performing linear canonical transformations in 
double phase space, we are free to choose the alternative coordinate planes, Q = (g_, q + ), 
or Q = x, or Q = y = J£ among others. In all cases, it is the fact that 

«W» = S(Q' ~ Q) (7-1) 
which permits us to identify the expansion coefficient in 

A=\A)) = JdQA(Q)\Q}} (7.2) 

with ((Q\A)). 

Let us now assume that the phase space is a product of a pair of phase spaces, 
X = X 1 £g> X 2 , each with 2Lj dimensions, and that these correspond to Hilbert spaces, 
Hj, so that H = Hi ® Hi- Then we can always decompose the Lagrangian planes 
chosen as a basis for double phase space as the product Q = Qi ® Q 2 , corresponding to, 
operators \Q)) = \Qi)) <8> IQ2)). Thus the complete \Q)) representation becomes 

A = \A)) = J dQxdQ 2 A(Qx, Q 2 ) \Q X )) <g> \Q 2 )). (7.3) 

The definition of the partial trace is then 

tr 2 A = tr 2 T 2 A = J dQxdQ 2 A(Q 1; Q 2 ) IQ^) «I 2 |Q 2 », (7.4) 

so that 

Ai(Qi) = J dQ 2 A(Q 1 ,Q 2 )((I 2 \Q 2 }} (7.5) 

defines the \Qi)) representation of a reduced operator A\, which acts on the Hilbert space 
Hi. It is well known that in the case of the density operator, p, the reduced operator, 
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pi, describes the same probability as the full density operator for all measurements 
concerning the subsystem- 1 . 

The different forms of the partial trace depend essentially on the Hilbert-Schmidt 
product (15.11) of each basis with the identity. In the case of the position basis, we have 

«I|Q)> = tr I |q_)(q + | = <5(q_ - q + ), (7.6) 

so that 

Ai(Qi) = J dg 2 _dg 2+ A(Q X , Q 2 = q 2+ )) <5(g 2 _ - q 2+ ) = J dq 2 A{Q 1 , (q 2 , q 2 )). (7.7) 
Here, we should recall that, 

A(Qi, (<?2, q2}) = (qi-, qi\A \ q 1+ , q 2 ), (7.8) 

in matrix notation. 

In the centre representation, we have simply 

«/|x»=tr/ (2 L R X ) = 1, (7.9) 

leading to the phase space projection: 

Ai(xi) = J dx 2 A( Xl ,x 2 ). (7.10) 

In the case of the density operator, the corresponding reduced Wigner function, Wi(xi), 
is thus obtained from W(x) in the same way as a marginal probability distribution is 
projected out of the full distribution (I3.12p . 

The simplest choice turns out to be the chord representation. Then, |/)) = T^ =0 is 
an element of the operator basis, so that 

«/|y» = S(y) = 5(0- (7.11) 

Thus in this case, instead of projecting, we obtain the reduced operator merely by slicing 
through the chord symbol: 

A 1 (6)=A(6,6 = 0). (7.12) 

Of course, the reduced operator A 1 itself is insensitive to the procedure used to obtain it 
within the various representations, but the ease of calculating the reduction is a special 
bonus of the chord representation. 

It should be recalled that the partial trace is invariant with respect to unitary 
transformations performed internaly within the factor Hilbert space H 2 : U = U 2 ® l\ 
(see e.g. [6]). In the example where the subsystems are particles that have separated 
by a large distance, then these are truly local transformations. In other words, if 
A' = UAlI -1 , then tr 2 A' = tr 2 A. This invariance corresponds semiclassically to the 
freedom of performing canonical transformations which leave invariant the xi variables: 
(xi,x 2 ) — > (xi,x 2 ). This also implies that only the double phase space corresponding 
to x 2 changes: (Xi,X 2 ) — > (X\,X' 2 ). If the canonical transformation is linear in the 

* A measurement on subsystem-2 only affects pi if the information on the outcome of the measurement 
is made available [5] . 
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single phase space, then both the centres, x, and the chords, £, are propagated in the 
arguments of their respective functions by this same transformation. 

Another point that is worth discussing concerns the completeness of the operator 
representations. Notice that the restricted translation operators 

| yi )) 1 = 4 =4<g>/ 2 (7.13) 

are a subset of the translation operators used in the chord basis for the full Hilbert 
space, Hi ® H 2 . It follows that a representation in terms of the restricted translation 
operators, , would not be complete. Likewise, we may define the restricted unitary 
reflection operators, 

|x 1 )) 1 = 2 Ll i?; i =i? Xl ®J 2 , (7.14) 

but these do not belong to the centre basis for H\ ® Hi- Even so, we may also define 
directly the reduced operator A 1 as 



Ax = J dxi Ai(xi) |xi))!, (7.15) 

with 

A 1 (x 1 )=trA(2 il ^ i ). (7.16) 

Let us now specialize to density operators. In the case of the chord function, we 
must take care of the normalization, which depends on the number of degrees of freedom. 
Hence, the validity of (17.171) between the chord representation of density operators, p(£) 
and pi(£i), implies that the reduced chord function is 

Xi(6) = (27rfi) L2 x(6,6 = 0). (7.17) 

Clearly, Xi(Ci) is the Fourier transform of Wi(xi). Since the definition of phase space 
correlations ( 16. 16ft is valid for the reduced system, we obtain the reduced correlations 
as a projection of the correlations of the entangled pure state: 

Ci(fi) = tr pi % j0) fa fj 1>0) 

= (2nH) L ^ J d Vl j^ h \ X i(Vi)\ 2 

= (27rfi) Ll J drnd V2 e^ A ^ h 5( V2 ) \ X (r])\ 2 

, d ^ f C(0 . (7.18) 
(2vrfi) L 2 KS1 K ' 

It should be recalled that the relation between the Wigner function and the chord 
function mimics that between a classical probability distribution and its characteristic 
function. The definition of correlations and the classical marginal distributions also 
goes through as above. Therefore, the property that the correlation of the reduced 
state for a given displacement, £i, is just the integral over all correlations in the larger 
space over displacements that share this component also holds for classical probability 
distributions. This relation does not depend on the full density operator being a pure 
state. 
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All the representations that we have been discussing will factor in the case that 
P = (l^i) ® | "02)) (("021 ® ("01 1) is a product pure state. Thus we obtain product Wigner 
functions, W(x) = Wi(xi) W^^) and product chord functions, x(0 — X2(^)- 
These relations may be interpreted in terms of average values of the basis operators, i.e. 
(Rx) = (-^xi) (^x 2 ) an< ^ (^) = (-^i) C^)- Thus, a sufficient criterion for the existence 
of entanglement would be that either of these equalities not hold for some centre, x, or 
some chord, £. 

Curiously, it is not the generation of cross correlations that is usually taken as a 
measure of entanglement, but instead the loss of correlations of the reduced density 
operator. Its von Neumann entropy (16.231) is often referred to as the entanglement. 
Expanding this to first order, results in the linear entropy of of a partial trace of the 
full density operator, 

l-trpi = l-d(0), (7-19) 

recalling (16.161) . which is the square of the concurrence, another widely used 
entanglement measure (see e. g. [59]). This is not an obvious measure of overall 
entanglement, because we should obtain the same measure by singling out instead the 
reduced density operator for subsystem-2. But, it is a simple consequence of (I6.24p . the 
invariance of the quantum correlations with respect to Fourier transforms for a pure 
state, that 



(2vrfi) J 

d6 f dv C( V )e^ 



(27Tfi) L2 J {2ttH) l 
d?7i 



(2 7 r^^ 1 '° ) = C72(0) - (7 - 20) 
Reinterpreted in terms of Wigner functions, 

J dx a [^(xO] 2 = J dx 2 [W^ 2 (x 2 )] 2 , (7.21) 

this is another remarkable property of pure quantum states, for it is highly unusual for 
the second moment of a pair of marginal probability distributions to display a similar 
equality. Indeed, it is not even generally true for product distributions. 

The focus on properties of the reduced density matrix makes sense when it is 
recalled that the concept of entanglement involves separate measurement on each of the 
components. The invariance of the partial traces with respect to local transformations 
carries over to the above measures of entanglement. Even better, it has been shown 
that it is possible to concentrate the entanglement within a few elements of an ensemble 
of identical states, by performing local measurements [60J. 

In terms of Husimi functions (16.381) . it is natural to describe entanglement in terms 
of the Wehrl entropy (16.391) for the reduced density operator. Another way of describing 
entanglement is through the Schmidt decomposition ( 13.18D . The corresponding Wigner 
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and chord functions are then 

W(x) = J>Aj (7-22) 

and 

X (0 = {2^)- L Y.K\ ^\\TL^\){^\\f'_^ 2 ), (7.23) 

recalling the definitions of the restricted reflection operators (17.141) and the restricted 
translation operators (17.131) . In both cases the partial trace over subsystem-2 substitutes 
the second Dirac bracket by Sij, so that 

W^xO = (vr^-^^A, 2 (^il^J^i) = (7.24) 

i i 

and 

XiCfi) = (27rfi)-^EA, 2 (^il^J^i) =EV^(6)- (7.25) 

Therefore, the reduced density operator is just a mixture of the factor states in the 
Schmidt decomposition for subsystem- 1, with probabilities specified by the square of 
the Schmidt coefficients. The concurrence is then given by 

l-tr(fc) 3 = l-5>*, (7.26) 

i 

in terms of the second moment of the weighing factors for the mixed state. Note 
that, contrary to the Schmidt number, this is a well defined entanglement measure 
for systems with infinite Hilbert spaces, if the above sum converges. Clearly, the purity 
of subsystem-2 involves the same sum over Schmidt coefficients, in agreement with our 
previous calculation (I7.20p . 

Consider now the case that a subsystem can again be split up into a pair of 
components. If the full original state was entangled, the reduced density operator is not 
pure. Hence, it is an average over pure states. Obviously, this cannot be a product state 
overall, but if all of the pure states are products, the mixed state is not characterized as 
entangled, rather it is a separable state. The problem with mixed states is that the 
decomposition into pure states is not unique, so a state is considered separable if there 
exists any decomposition where it is separated. 

Let us now define a classical pure state as a 5-function in phase space. Then all 
pure states, f(x) = S(x), in a higher dimensional phase space will be product states, 
because the higher dimensional (^-functions factor as S(x) = S(xi) <5(x 2 ). In this sense, 
the expression 

f( x ) = J f( x ') 5(x - x')dx' (7.27) 

can be reinterpreted as a classical separable state: Any probability distribution in phase 
space can be considered as a linear combination of products of classical pure states. 
Thus, we can never consider a classical phase space distribution to be entangled, no 
matter how strong the correlations may be between variables pertaining to different 
subsystems. 
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What if a mixed Wigner function for a bipartite state is everywhere positive? Can 
we mimmick the above reasoning to conclude that there is no entanglement? In general 
this is not so, because the function S(x) does not represent a density operator in the 
Weyl representation. It represents instead the reflection operator, which has an infinitely 
degenerate negative eigenvalue, as discussed in section 5. The closest that is possible is 
the coherent state (16.301) . which approaches a 5-function as H — > 0, but imposes an extra 
smoothing on the Wigner function for any combination of these pure states. Indeed, a 
general superposition of coherent states is defined by a weight function, known as the 
Glauber-Sudarshan P-function in quantum optics [611 [621 [8] . So, it is the positivity of 
a P-function that guarantees a separable state, rather than that of the Wigner function, 
because each coherent state can be factored. 

To close this section, let us now study another kind of projection of the Wigner 
function. Whereas, by projecting onto a component subspace we generate a reduced 
Wigner function, a projection onto a Lagrangian plane (16. 9p results in a probability 
density. All the coordinates of such a plane correspond to commuting operators. In the 
case of a bipartite system, we can define this Lagrangian plane by choosing some linear 
combination of the variables for each subsystem, q[ = a±pi + fiiqi and q' 2 = a 2 p 2 + /?2<?2, 
so that each coordinate, q'j, pertains to a diferent subsystem. 

Consider now pairs of either-or measurements on both these variables, i. e. we can 
define observables Oi a , 0\ b , 2a and 2b which take the value +1, for q'j in the interval 
ja, and —1 outside. In terms of projection operators Pj a , we have Oj a = 2Pj a — I and 
OxaOia = ^P\aP-ia — 2Pia — ^P 2a + I, with similar formulae for the other products of 
commuting operators. Combining the expectation values for these products in the form 
of the CHSH inequality fEUTL 

(d la d 2a ) + (d la o 2b ) + (d lb d 2a ) - (d lb d 2b ) = A\(P la P 2a ) + (P la P 2b ) + (P lb P 2a ) - (Pi b P 2b )} 

- 4[(P la ) + (P 2a )\ + 2, (7.28) 

we can now evaluate each expectation value on the right hand side as a definite integral 
of the probability density over some region of the (q[, q' 2 ) plane. This is a purely classical 
setup, so that by regrouping, 

2 - 4[P la - (P la P 2b )] - 4[P 2a - (P lb P 2a ) - (P la P 2a )} - A(P lb P 2b ) < 2, (7.29) 

we rederive the CHSH inequality, because the square brackets cannot be negative. 

We thus verify that the correlations measured among commuting pairs of either-or 
observables of each subsystem lie within strictly classical bounds, irrespective of any 
possible entanglement of their combined state. It makes no difference whether, or not, 
the Wigner function has negative regions. The point is that we need only deal with a 
single positive projection, which is a true probability distribution. To obtain a violation 
of the CHSH inequality, we must choose pairs of observables for each component which 
do not commute. The correlation for a given choice of observables from each pair may 
still be computed from the probabilities in the respective Lagrangian plane, but we must 
use different planes in each of the four correlations. Then, if the overall Wigner function 
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that generates all these densities has negative regions, the CHSH inequality may be 
violated, as discussed in the following section. 

Apparently, there has not been much effort to relate the intuitively appealing 
picture of entanglement as the source of nonclassical correlations in Bell inequalities 
to the technical entanglement measures appropriate to quantum information theory. 
However, a recent paper by Cirone [63] bridges this gap for finite dimensional systems. 
The main point is that measurements are restricted to projection operators for the 
factor states in the Schmidt basis. It is then shown that the same concurrence, which 
was introduced in terms of the partial trace, can be expressed as a sum over conditional 
probabilities for measurements on either component. 

8. Generating a classical entanglement: The EPR state 

We have seen how symplectic transformations correspond exactly to unitary 
transformations in Hilbert space. Let us now examine how these can produce entangled 
states, given that the initial state, p, is a product of states, each represented by its 
Wigner function, W^Xy), or its chord function, Xj(£j)-> so that W(x) = W\ (xi) W 2 ( x 2) 
and x{£) = Xi(^i)X2{^2) ■ For the canonical transformation to be linear, the classical 
interaction Hamiltonian, H{x\,Xq) can only be bilinear in the phase space variables. A 
convenient choice is H = piq 2 — P2Q1, which may be interpreted as angular momentum, 
L3, if the degrees of freedom refer to Cartesian coordinates in a plane. This Hamiltonian 
merely rotates both p and q coordinates in the argument of W(x) and x(Q- Then, after 
a rotation by 7r/4, the density operator becomes p / , represented by 



Since the partial trace is specified by (17.171) . a section of the chord function, the reduced 
density for the first component becomes 



in the chord representation. 

So as to emphasise how classical an entanglement can be, let us choose for example 
initial Gaussian states, the product of harmonic oscillator ground states, described by 




(8.1) 






1 



1 



(8.3) 



or 




(8.4) 



Thus, the probability distribution for positions, 
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is also a Gaussian with elliptic level curves, that are also rotated if u)\ 7^ uo 2 . In this case, 
the effect of rotation, followed by the partial trace, is just a narrowing of the Gaussians 
in the chord representation: 

2 , ~ - . />. \ 2" 



xi(6 



2nh 



exp 



2h 





(8.6) 



Notice that normalization is maintained, because we still have Xi(Ci) = (2iTh)^ 1 at 
the chord origin, but now the widths of the position Gaussian and of the momentum 
Gaussian are obtained through different kinds of average. The overall narrowing 
indicates that this is no longer a pure state. 

The Wigner function presents a more intuitive picture of a mixed state. Taking the 
Fourier transform: 



W i(*i) = T ex P 

7TA 



2 P ? 



(8.7) 



H(uji + tu 2 ) h(ui + u 2 ) 

This still integrates to one, as demanded by normalization, but the Gaussian is now 
broader, with the uncertainty A = {yji + u 2 ) /2^/ujiuj 2 > h, if oj\ 7^ u 2 . Therefore, this 
is not a pure state. The way that this example relates entanglement to initial states 
and evolution, which may both be considered classical, is even more extreme than those 
discussed in [67], which relie on projections of the Husimi function, in the approximate 
role of phase space probability density. 

Another confirmation that this is not a pure state is that 

tr (p^) 2 = 2nh J dxx [W^)] 2 = 2nH J d& |x'i(6)| 2 < 1, (8.8) 

and yet it might seem that this is just a freak result. After all, our state has remained a 
smooth classical-like Gaussian throughout. There are none of the quantum oscillations 
which are supposed to be the fingerprint of nonclassicallity: For a start, nothing would 
prevent us from identifying the original Wigner function with a classical probability 
distribution. We then perform a simple rotation with perfect classical correspondence 
and obtain a new Gaussian, which pretends to be a quintessentially quantum entangled 
state! Have we been led astray? 

Let us go back to the full Wigner function, resulting from the choice (18.31) of 
Gaussians for the initial product state. After the 7r/4 rotation, this is just 



W't 



-V 



ixh) 



exp 



~2H 



+ Q2Y ~ 



In the extreme limit where U\ 
of 



2hjj\ 
and uj 2 - 



(P1+P2)' 



cxp 



^2 , 



q-2, 



2hioo 



(Pi ~ Pa) 



.9) 



00, we would obtain a normalized version 



W'{x)=6(q 1 -q 2 )6(p 1 +p 2 ), (8.10) 

which is just the Wigner function derived by Bell [Mj for the original EPR wave function 
[3], namely (q\if)) = <5(<?i — 92) ■ It so happens that the rotation that transformed the 
coordinates of our initial state, i. e. the ground state of an anisotropic plane harmonic 
oscilator, is essentially the same as the transformation from the individual coordinates 
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for a pair of particles into a centre of mass, together with a relative internal coordinate. 
(These transformations differ only by local unitary transformations.) The EPR state 
is a ^-function both in the relative position and in the total momentum, which is the 
conjugate variable to the centre of mass. 

Thus, the entanglement verified in our initial example implies that the centre of 
mass is likewise entangled with the relative coordinate in the EPR state. Perhaps, it is 
then even more surprising that the example that was picked is in some sense classical, 
if we consider that the discussion of the nonlocal and hence nonclassical nature of 
quantum mechanics started off with the historic EPR paper [3j. The fact that the full 
Wigner function is positive, not only allows us to interpret it as a classical probability 
distribution, but it also ensures that there is a wide range of measurements that can 
be performed on either component which may be considered as classical and hence 
local. We already found in the previous section that any measurement of pairs of 
either-or variables, Oi a , Ou, 2a and 2 b which take the value +1, for general phase 
space coordinates, q'j, in the interval ja, and —1 outside, have correlations that satisfy 
the CHSH inequality, even if the Wigner function has negative regions. That was the 
case where the quantum observables which are measured commute. The statement 
for positive Wigner functions, due to Bell [64j, is stronger: The inequality is then 
satisfied even if we choose different variables for each measurement, q' la ^ q' lb and 
q' 2a 7^ q' 2b , corresponding to different Lagrangian planes in phase space and, hence, 
quantum operators that do not commute. The argument is essentially the same as in 
the last section, except that now we can obtain all the expectation values from the 
full Wigner function, acting as a global probability distribution, instead of dealing with 
different probability distributions, each restricted to the Lagrangian plane specific to a 
given pair of variables. 

Let us now reexamine our classically entangled states from the point of view of 
the reduced reflection operators, R'^j, defined as (17.161) . that act on each component 
and, in particular, the parity operators, RL. Such observables do not correspond to 
smooth phase space fuctions in classical mechanics, indeed, the Weyl representation of 
these operators (15.301) is singular. Nonetheless, parity, or reflection measurements can 
also be carried out on classical waves, as discussed in section 2, and the question now 
concerns the possible correlations between measurements for different reflections carried 
out on both components. The fact that the full Wigner function (I8.9P is symmetric 
with respect to the origin implies that the density operator commutes with the full 
reflection operator, R . However, W((0) < nh, so it does not have pure parity, i. e. p 1 
does not commute with R' Ql . Hence, according to the discussion in section 6, there is a 
finite probability to obtain negative (odd) parity, if such a measurement is performed 
on subsystem- 1. 

Perhaps this would not be so obvious a priori: The original state, represented by 
Wo(x), is a pure state with pure positive (even) parity and this is also a property of 
the rotated state. This property can be verified directly, or it may be noticed that the 
driving Hamiltonian commutes with R , so that H(x) = if(i? (x)). But now we find 
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that a measurement of the parity of subsystem-1 has a finite probability to be negative. 
How is that? 

Notice that the same also hods for subsystem-2: The derivation of the reduced 
density operator, p^, represented by W^^) an d Xzi^), goes through exactly as above. 
Therefore there is also a finite probability of measuring negative parity in subsystem-2. 
As was shown in section 6, the fact that, in both cases, the Wigner function is symmetric 
about the origin implies that all the pure states, into which the mixed reduced density 
operator can be decomposed, must have pure parity, but they are not all even. For this 
reason, the Wigner function (18. 7p had to be obtained as a Fourier transform of the chord 
function; not a mere rescaling. 

The crucial point is that the rotated state, p 7 , does not commute with either of 
the restricted reflections defined by ( 17.141) . i.e. R' Q1 or R' 02 , even though it commutes 
with their product: Rq = R' 01 R f 02 = R^R'oi- ft should be recalled that the reduced 
Wigner functions are entirely determined by (17.151) and (17.161) in terms of the restricted 
reflections. Thus, to understand the results of measurements of either R' j, we need a 
common basis for all these operators. This is just the product of an even-odd basis for 
subsystem-1 and subsystem-2, for which we obtain the table: 



even ® even — > even (8-H) 

even <8> odd — > odd (8-12) 

odd <S> even — > odd (8.13) 

odd (g> odd — > even (8-14) 



Since fl is even, it must be a superposition of the subset of basis states: even eg) even, or 
odd® odd. Furthermore, we now find that the evolved state has a full parity correlation: 
If the measurement of R' Q1 specifies even parity, then this must be the outcome of a 
measurement on R' 02 . Conversely, if one of the subsystems has odd parity, then we 
know this to be the parity of the other subsystem. 

An initial product state of an even Schrodinger cat state with a coherent state, 
which is rotated by tt/4, is also susceptible to the foregoing analysis. However, an odd 
symmetry Schrodinger cat would have perfectly anticorrelated odd-even, or even-odd 
subsystems. In the case of the rotated cat the evidence for entanglement is much more 
obvious. The pair of Gaussians is not centred on either of the planes in the chord phase 
space pertaining to the pair of subsystems. The partial trace that generates the reduced 
chord functions is a section of the full chord function, so that it does not capture these 
local maxima. Therefore, there is a deficit of phase space correlations in the reduced 
density operators. 

Returning to the original rotated squeezed state, or, equivalently, the original EPR 
state, we must conclude that this is truly quantum and correctly described as entangled, 
i. e. just as nonclassical as the spin states in the Bohm version of EPR [7J that are 
commonly used to exemplify entanglement. The secret lies in choosing the property to 
be measured: A position measurement on one of the subsystems would not distinguish 
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between this pure quantum state and a classical distribution. However, a measurement 
of reflection eigenvalues evokes a spin-like duality of this apparently classical state. 

The violation of the CHSH inequality for reflection measurements of the smoothed 
EPR state completes the evidence of its nonclassicality. Banaszek and Wodkiewicz 
[65] first pointed out that the full pure state Wigner function of a bipartite state 
is proportional to the correlation for relection measurements on each subsystem: 
(R'xiR'ia) = (7r7i) 2 W(x). This leads to a violation of the CHSH inequality for reflection 
measurements of the EPR state. They also proposed a realistic experiment for this 
in quantum optics |66j. We have already verified the complete correlation for parity 
measurements about the origin, which is in agreement with the maximal value that 
the full Wigner function (18.91) attains there. Its decay for large X! or x 2 signifies that 
(K1K2) 0, so that 

Cchsh = (R' 01 R' 02 } + (#01^x2) + (RxiR'02} - {R'*iR'x2) (8.15) 

sinks from, 2, its maximal classical value at the origin to the limiting value 1. However, 
the origin is not the maximum of Cchsh^i, x 2), because the lowest order expansion of 

(7rfi) 2 H/'(x) = 1-^ (Vi + K) 2 -^- (Pi+P2) 2 -Yn ^~^ 2 -^tr 2 (Pi-P 2 ) 2 + -(8.16) 
leads to 

Cchsh{*i, x 2 ) = 2 + (uji - uj 2 )qiq2 + (— -) P1P2 + ••• (8.17) 

\UJl L0 2 / 

Hence, the origin is a saddle point of Cchsh^i, x 2)> which increases from its maximal 
classical value along the directions q 2 = —qi and p 2 = pi, if one chooses the EPR 
conditions, uj 2 » u\, i.e if reflections are chosen in the directions where the Wigner 
function decays rapidly. 

So we find that nonlocal correlations between two subsystems can arise even if 
the Wigner function for the full system is everywhere non-negative. It would thus 
appear that there is no relation between fringes in the Wigner function, where it 
attains negative values, and entanglement. The former project as interference fringes 
for possible measurements, but this is quite a different kind of nonclassicality than 
the delicate nonclassical correlations resulting from entanglement. But even here, one 
must be wary! If the measurements on the different components concern mechanical 
observables, natural for classical particles, then there is at least one case where negativity 
of the Wigner function has been shown to produce nonclassical correlations. Indeed, Bell 
[64"] constructed an example where the CHSH inequality is violated for measurements 
on pairs of different variables, q'j = qj + tjpj and q'j = qj + Tjpj. The state for which 
this is proved is a variation of our rotated state, where one of the factor Gaussians is 
substituted by the second excited state of the harmonic oscillator. 

It should always be remembered that entanglement is not an intrinsic property, 
but only acquires its meaning within a specified basis, the computational basis, or 

jj Even though Bell reffers to the transformation parameters as times, these should be understood as 
specifications of the variables and hence of the planes onto which the Wigner function is projected. 
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the basis where measurements are made. In this respect, it resembles semiclassical 
caustics, which depend on our choice of representation. If the physical realization of the 
foregoing example were the ground state of a 2-D harmonic oscillator, then the rotation, 
which was found to produce entanglement, could be dismissed as merely an inconvenient 
coordinate transformation: Unless all measurements were to be restricted to the original 
coordinate axes, it would not be relevant, though true, to say that the rotated system 
became entangled, while the original system was a mere product. In contrast, for the 
alternative physical interpretation of one of the new coordinates as the centre of mass 
for a pair of particles, its entanglement with the internal coordinate can be important. 

9. Entanglement and decoherence 

The process of decoherence also results from the interaction of a pair of systems; the 
(small) open system and a (large) system, which we call the environment. In contrast to 
the previous example, the component over which we trace, so as to obtain the reduced 
density operator, is on a scale which defies anything but a statistical description. The 
ususal picture is that the environment lies somewhere outside, but it may just as well 
consist of the internal degrees of freedom for the centre of mass (CM) of a large system 
of particles. Exchanges between the large scale motion and the internal variables lead 
to macroscopic energy dissipation as well as decoherence of the quantum state for the 
CM. 

Let us consider the simplest possible example of the decoherence of the CM, because 
of its entanglement with internal variables. The CM for a system of L identical particles, 
assumed to be distinguishable is Q = (qi + ... + q£)jL. The conjugate variable to Q is 
the total momentum, P — pi + ...+pl. Let us further imagine that they are each in the 
same single particle state, p, and that these are independent, i. e. both the Wigner and 
the chord function are products over those of the individual states. This may seem too 
restrictive, because we should allow for different values of each average position (qj), 
but we can redefine this as the origin for each j, so that we then measure Q from (Q). 

In the case of L = 2, X = (P, Q) is obtained from the rotated coordinate in 
the previous section by a mere canonical rescaling of 2 ±1//2 . It has been repeatedly 
emphasised that all such symplectic transformations on the argument of the Wigner or 
the chord function correspond exactly to unitary quantum transformations. So let us 
now reverse this transformation in the case of general L: We define Q' = L l l 2 Q and P' = 
L V I 2 P. Then, if the individual Wigner functions, W(xj-), were classical probability 
distributions, the Central Limit Theorem would imply that the distribution for 
X' = (P', Q') converges to 

M^(X') -> [TrAKrW-X'K^Xy^, (9.1) 

as L — > oo, where we recall the definition of the Schrodinger covariance matrix, K, 
in (13.71) and its determinant Ak 2 . It is remarkable that positivity is not a necessary 
ingredient for the proof of the Central Limit Theorem: In the case of identical square- 
integrable pure state Wigner functions, it is shown by Tegmark and Shapiro [70], that 
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convergence onto a Gaussian again results. If the state for the individual particles is not 
represented by a pure state Gaussian, then the moments for this state will be such that 
Ak > h. Therefore, the centre of mass Wigner function, Wl(X), is a broader Gaussian 
than is permissible for a pure state and hence it must be a mixture. So, the CM of 
independent particles with identical Wigner functions is generally entangled with the 
internal phase space variables (which it has not been necessary to describe explicitly). 
Curiously, the potential entanglement resulting from the Central Limit Theorem was 
overlooked in [70] . 

How does this entanglement with the internal coordinates evolve in time? It is easy 
to verify that free motion, generated by the Hamiltonian H = pi 2 + ...+pL 2 , will not alter 
Ak- Thus the entanglement of the centre of mass with the environment is invariant in 
this simple case. Let us supose instead that, though the particles do not interact, there 
is an external nonlinear field. Furthermore, the particles are sufficiently separated and 
the field is smooth enough, so that it is legitimate to linearize the field locally around 
each (qj). Then the Wigner function for each particle will evolve classically in different 
ways. The restrictive form of the Central Limit Theorem in [70J cannot be applied 
in this case, but one can readily adapt Levy's proof [71] to allow for different Wigner 
functions, as long as the moments are finite and their average values congerge [72J. 

The averages of the moments resulting from the different evolutions of many Wigner 
function lead to a progressive loss of purity for the CM. Just as in (18. 7ft for the simple 
example of the last section, the uncertainty, Ak, increases. On top of that, the Central 
Limit Theorem supplies the statistical ingredient for the decoherence process. It might 
appear strange to obtain decoherence even for a system of noninteracting particles, but 
it should be recalled that the CM momentum P, or P', appears linearly in each of 
the terms, Pj 2 , in the kinetic energy, which accounts for the coupling to the internal 
momenta. 

So as to make contact with the theory of Markovian open systems, we can now 
reinterpret this evolution of the reduced density matrix as a convolution of the original 
(Gaussian) Wigner function for the CM with a broadening Gaussian. For its Fourier 
transform, the chord function, this evolution is merely the product of an initial Gaussian 
with another Gaussian that narrows in time. This is exactly the result for quantum 
Markovian evolution of an open system, in the case of quadratic internal Hamiltonian 
and linear coupling to the environment |73j . 

The deduction of the canonical Lindblad equation (see e.g. |74j ) . 



that governs the evolution of the density operator in the quantum Markovian theory does 
not proceed by tracing out a larger system. All the same, the mere fact that the evolution 
is entirely determined by a differential equation precludes any delayed participation of 

ff It must be recalled that the average CM evolution, (X(t)) has been hidden by a time dependent 
coordinate transformation. 



dp 

di 



(9.2) 
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previous motion. The Lindblad operators, Lj, account for the nonunitarity of the 
evolution, that is, they take the part of the coupling to the environment. The Markovian 
approximation can in principle include arbitrary (non-quadratic) internal Hamiltonians 
for the system. 

The derivation of the Markovian approximation in the context of quantum optics 
(the damped harmonic oscillator) was carried out originally by Agarwal [76J , but this is 
all in the language of complex phase space. The exact solution of (19.21) in [[76] and that of 
Diosy and Kiefer [75], for the free open particle, are special cases of of the general result 
in [73]: The chord representation of (19. 2\i is particularly simple if the Lindblad operators, 
are linear functions of positions and momenta, Lj = lj ■ x, and if the Hamiltonian is 
quadratic [73] : 

t) = {H®, x& *)}"^E & ■ 2 Xfr t). (9.3) 

3 

Here, the first term is the classical Poisson bracket. The exact solution of this equation 
factors into the unitary evolution of the chord function, undistrubed by the Lindblad 
operators, and a narrowing Gaussian factor. In the Wigner representation this becomes 
a Gaussian smudging of the unitarily evolving Wigner function. It is remarkable that 
the Wigner function becomes positive after a time that depends only on the parameters 
of the Lindblad equation, regardless of the initial pure state (75j [73] . 

In our simple example of the evolution of the CM, the Lindblad operator for its 
one-dimensional motion, should be chosen as the total momentum P, because this is 
the variable that couples to the internal motion, which is hidden within the Markovian 
approximation. Even though the Central Limit Theorem supplied the Gaussian factor 
of the evolving chord function, the overall Gaussian form for the evolving CM does 
not reflect the richness of other possibilities for Markovian evolution. However, by 
considering the entanglement of a small system with the CM of a large system and 
following the treatment of the example in the preceding section, we obtain qualitatively 
the general Markovian picture. This allows an interpretation of the Gaussian smoothing 
as originating in the multiple small contributions contemplated in the Central Limit 
Theorem. 

The standard way of going beyond the Markovian approximation, so as to include 
memory kept by the environment of the previous motion of the system, is to use the 
Feynman- Vernon functional [77], in the manner exploited by Caldeira and Leggett [78J. 

10. A semiclassical picture of entanglement 

A full semiclassical theory of entanglement is still a program for the future, fascinating 
but difficult. However, several of the main elements are sketched in this concluding 
section. 

For a start, one should note that it is feasible to fit semiclassical torus states with 
Gaussian coherent states placed along the classical torus in a very satisfactory way 
[79] . The number of Gaussians required increases with a fractional power of The 
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important qualitative feature is that the interference fringes of the Wigner function, 
near the midpoint of the pair of Gaussians composing a Schrodinger cat, have the 
same wavevector as the similar fringes at the centre of a geometrical chord of the 
classical torus. Therefore, in both cases we can describe very fine interference fringes 
related to long chords. It also follows that our preliminary study of entanglement 
and decoherence of cat states is not at all irrelevant for understanding the evolution 
of product semiclassical states. Refinements of the fitting procedure allow even the 
description of the diffraction effects near caustics [80J. 

Before analysing product states and their partial trace, recall that dyadic operators, 
live in a kind of squared Hilbert space, which corresponds to a double phase 
space. These operators were shown in section 5 to correspond to a product Lagrangian 
surface in double phase space, r^,®r^, if each of these states corresponds to a Lagrangian 
surface on its own right. Thus, the projection operator, or pure state density operator, 
p = \if))(if)\, is just a particular case of this general rule. If the state, corresponds to 
a Bohr-Sommerfeld quantized torus of L dimensions in a 2L-D phase space, then the full 
density operator must correspond to a 2L-D product torus in 4L-D double phase space. 
This is in exact analogy to the way that a product torus describing the state for several 
particles (I4.12p is obtained from lower dimensional tori. Recalling that we can describe 
double phase space in terms of the centre coordinates, x, and the conjugate variables, 
y = J£ (15.1 2D . the semiclassical Wigner function, W(x), is then a superposition of 
complex exponentials, such that each phase is obtained by integrating y(x) along one of 
the different branches of the torus. Even though this approximation breaks down along 
caustics, the latter provide ready indication of regions where the Wigner function has a 
large intensity. 

The problem is then to relate the semiclassical Wigner function, defined on the 
centre plane, to classical structures that are also portrayed in this same single phase 
space. Let us consider first the semiclassical Wigner and chord functions in the simplest 
case where L = 1. The Fock states (I6.34p are good examples of semiclassical torus 
states when the quantum number n is large. Introducing the asymptotic expression for 
Laguerre polynomials, 



brings the Wigner function for these states into a semiclassical form. To understand 
this, we must investigate the geometry of the double torus from the point of view of the 
simpler quantized curve, which is just a circle in this case. 

Every point on the double torus represents a pair of points on the quantized curve 
and vice versa. A given pair of points on the quantized curve, x±, defines a geometric 
chord: £ = x + — x_. Hence, y = J£ is the chord coordinate on the double torus, which 




(10.1) 



together with the large argument expansion, 




(10.2) 
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Figure 8. (a) The semiclassical Wigner function is constructed by reflecting the torus 
r around the point x: The chords £(x) are defined by the intersections with R x t. (b) 
The semiclassical chord function is constructed by translating r by the vector £: The 
centres x(£) lie halfway back along £ from each intersection of t with T^t. 

has the centre coordinate, x = (x + + x_)/2. Obviously, the exchange of x + with x_ 
produces a new chord of the quantized curve with the same centre, x. Viewed in double 
phase space, there must always be pairs of chords of the double torus projecting onto 
each centre, x. The symmetry of this surface with respect to the identity plane, y = 0, 
leads to complex conjugate phase contributions, in line with the above cosine for the 
Fock state. Actually this is a general feature: Because the Wigner function is real, the 
chord pairs will always produce semiclassical contributions adding up to cosines. 

To obtain the phase of the cosine contribution to the semiclassical Wigner function 
for each pair of chords, the best course is to use a result which was put in its most 
general form by Littlejohn [31]. This concerns the general overlap, (ip\<f)), of quantum 
states associated semiclassically to curves and r^: The semiclassical contributions 
arise from the intersections of these classical curves and the phase difference between 
a pair of contributions is just the area sandwitched between the corresponding pair of 
intersections, divided by Planck's constant, as shown in FigJHJ 

We can immediately apply this principle to the Wigner and chord functions for pure 
states, by recalling that W^(x) = (tt^^K^xI^)) and xv(0 = ( 27rfi ) _1 MCW>)- 
The semiclassical state i2 x |^>) is merely the state constructed from the reflected curve, 
R*.( T ip), whereas T^\ip) corresponds to the translated curve T^{r^). Therefore, in the case 
of the Wigner function, we obtain the phase of the cosine as half of the area sandwitched 
between the torus and its reflection at the centre x [22], which coincides with the area 
between the torus and the chord [HUES] ( see a l so - [H3]). Furthermore, this construction 
supplies, at a glance, the tips of all the chords centred on x as the intersections between 
both curves, as seen in Fig.[S(a). Note that, once the curve has been reflected around the 
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origin, we need only translate Rq(t^) to obtain the reflections around all other centres, 
because of the group property, i? x = T x o R o T_ x . 

The same geometrical method can be used to study the structure of centres for a 
pre-specified chord, £, on the curve, r^. Each intersection of t$ with the translated curve, 
T^T^p), reveals one of the tips where the chord is to be placed and hence the centre of the 
chord, as shown in Fig.^b). In the case of open curves, the chord function may actually 
be simpler than the Wigner function, because it is not necessary to have interference, 
as in the case of the parabola, for which there is only one intersection. In the case 
that T^p has a centre of symmetry, as in the example of the Fock state, we thus find 
that the simple relation (16. 29ft between Wigner and chord functions is respected by the 
semiclassical approximation. 

Viewed in single phase space, the caustics of the Wigner function arise from 
coalescing torus chords, as their centre, x, is moved. This occurs at the tangencies 
of -Rx( r i/)) with the fixed curve, r^. Similarly, the caustics of the chord function are the 
loci of £ such that T^(r^) is tangent to t$. On the other hand, in double phase space, 
the Wigner caustics for a torus state are viewed as projection singularities of the double 
torus, T^p <S> T^p, which lies above the area inside r^. The general geometric constructions 
underlying the semiclassical Wigner and chord functions are readily extended to phase 
space representations of dyadic operators, \ip){(fi\, corresponding to double tori, ® t$. 
Their Weyl representation is known as cross-Wigner functions or Moyal brackets 
[68] . whose semiclassical form is presented in [69J. 

There are many fascinating features of caustics in the phase space representations 
of pure states that have been studied and many more that must still be analyzed. For 
instance, the build up in the centre of the Wigner function for the Fock state is a caustic. 
Its semiclassical origin is the degeneracy of a continuum of chords conjugate to the same 
symmetry centre. However, this is a nongeneric feature of reflect ion- symmetric states. If 
the symmetry is broken, this supercaustic unfolds into a cusped triangle, first described 
by Berry [82]. The unfolding of higher dimensional caustics for rotated product tori, 
studied in [22] were also examined for the Wigner function. It turns out that the double 
fold surfaces of the Wigner caustic that meet along the torus do not unfold in the manner 
portrayed in Fig HI because of a symmetry constraint. 

The limit of small chords is specially relevant for semiclassical theory. For the 
Wigner function, it singles out the classical torus itself as the Wigner caustic. The 
uniform approximation for the Wigner function throughout this region is presented in 
[82] , for the case of a curve and in [84] , for a two dimensional torus. A pair of sheets 
of the double torus are joined on the identity plane along this curve, or torus. The 
large amplitude of the Wigner function oscillations near the quantized curve is due 
to this caustic. The corresponding caustic of the chord function collapses onto the 
origin, whatever the geometry of the classical region. The neighbourhood of this highly 
nongeneric chord-caustic is discussed in jl2]. Once again, we find that all regions where 
C(£) = |x(0| 2 i s large, outside of an fi L//2 -neighbourhood of the origin, point to phase 
space correlations that are truly quantum in nature. 
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Figure 9. Typical fringes for a scmiclassical Wigner function. The classical curve 
follows closely the border of the fringes. The triangular structure of interfering fringes 
near the centre results from the umfolding of the maximum of the Fock states due to 
symmetry breaking. 
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So far, we have only discussed static properties of the density operator. In turning 
to dynamics, a preliminary point is that we should distinguish the Weyl propagator, 
U(x), that is, the Weyl representation of the evolution operator, from the propagator 
for Wigner functions. The former is unitary and is hence supported by its own 
static Lagrangian surface in double phase space, as discussed in section 5, so that its 
semiclassical description is similar to that of the Wigner function itself. This can be 
deduced from a path integral in single phase space [85], [31] . 

So as to treat the unitary evolution of the density operator for a pure bound state, 
\ip)(ip\ ) with L degrees of freedom, we need to consider the corresponding classical 
evolution of a 2L-D Lagrangian torus, ® t$. Initially, this is separable within 
both single phase spaces, even though the product is not factored in the centre x chord 
coordinates. The classical motion must propagate the tips of each chord, x_ and x+, 
in the same way. Taking account of the change of sign, p_ — > — p_, in the definition of 
double phase space, we find that the double phase space Hamiltonian must be 

H(X) = H(x+) - H{xJ) = #(x - Jy/2) - H(x + Jy/2). (10.3) 

This classical Hamiltonian can be verified to preserve the product form of the 
geometric structures in each of the phase spaces x±, but it will not preserve initial 
products within each of these in the general case that the single Hamiltonian H(x) has 
coupling terms between different degrees of freedom. It propagates Lagrangian surfaces 
in double phase space that correspond either to density operators (according to the 
Liouville-von Neumann equation), or to unitary operators (the Heisenberg equation). 
The explicit formulae for the semiclassical evolution of the Wigner function are given 
in [86j EZ] , whereas the evolving action of the chord function is presented in [88] . The 
difficulty that cannot be avoided by changing the representation lies in the caustics of 
the initial state, which require more sophisticated semiclassical treatment. 

A promissing approach lies in the definition of integral propagators for the Wigner 
function, or for the chord function. The former may be defined in terms of the Weyl 
propagator as a kind of second order Wigner transform (see e. g. [86] ) 

U(x t ,x)= /d/i?7_ t (^^-^)C/ t (^^ + /i)exp(^//A(x t -x). (10.4) 
J 2 In 

Their explicit semiclassical form has been developed in [89], but these propagators 

also have their own intrinsic caustics. More recently, caustic-free propagators, from 

the Wigner to the chord function and vice versa, have been defined [88]. These are 

constructed either in terms of the propagation of the unitary reflection operators, or the 

translation operators, instead of directly evolving the density operator itself. 

Several of the geometrical structures underlying the semiclassical theory of the 

Wigner function for nonseparable tori in double phase space that evolve under the action 

of a general Hamiltonian were analyzed in [I22J for the simplest case where L = 2. It will 

be necessary to push much further this analysis, while adapting it to the chord function. 

The reason is again that the partial trace of the density operator is obtained immediately 

by a section through the chord function, £2 = 0), which is defined semiclassically by 
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a nonseparable section of the double torus that evolved from an initial product <g> r^. 
For instance, the quantum state and the corresponding double torus both loose their 
product form under the action of the simple Hamiltonian employed in section 8. But, 
slicing through a torus produces either a single, or several lower dimensional tori. 

This indicates that the semiclassical theory of reduced density operators that allow 
us to quantify entanglement, or to calculate the correlations on separate measurements 
effected on each component of the system, can still be associated to lower dimensional 
Lagrangian surfaces that are no longer products. The problem is to work out the 
actions and the amplitudes for this multidimensional geometry. This general picture 
agrees with initial results for nonunitary Markovian evolution of semiclassical Wigner 
functions (90]. It is notable that the same methods which have been used in [91J to show 
that the semiclassical approximation to Wigner function satisfy the purity condition, 
tr p 2 = p, reveal the loss of purity with time due to decoherence. 

A final comment about entanglent concerns the choice of components which arises 
with the freedom provided by unitary transformations, corresponding to classical 
canonical transformations. It is possible to disentangle each state that has lost its 
original product form by merely reversing time. This corresponds classically to runing 
Hamilton's equations backwards in the double phase space. In the example of section 8, 
the product of two harmonic oscillator ground states is recovered from the EPR state. 

In quantum mechanics, one can always specify a unitary transformation on the 
entire Hilbert space, which transforms any given state into any other state and we 
can choose the latter to be a product state. However, the correspondence for this 
disentanglement cannot exist for typical quantum eigenstates of a classically chaotic 
Hamiltonian. According to Shnirelman's theorem [92J, these ergodic states are supported 
by the full energy shell, in the sense that averages of smooth observables are well 
approximated by classical averages over this surface. In the case where L = 2, the energy 
shell is 3-D and, because it has an extra dimension, there exists no classical canonical 
transformation, whether linear or nonlinear, that can transform it into the product of two 
closed curves. Therefore, ergodic eigenstates are essentially entangled from the point of 
view of classical correspondence. As it happens, it is precisely this set of states that has 
resisted for decades all attempts at a compact semiclassical characterization. The study 
of traces of classical chaos in quantum mechanics is known as quantum chaology [83J. 
The characterization of ergodic states as those that are not classically disentangleable 
creates a bridge between entanglement theory and quantum chaology. It remains to 
be seen whether this special type of entangled state has any application in quantum 
information theory. 

Acknowledgments 

This text develops a sereis of Lectures delivered at the MPIPKS-Dresden in September 
2005. Many of the ideas and results presented here arose in the context of a fruitful 
collaboration with Olivier Brodier in CBPF. Raul Vallejos, Marcos Saraceno and Pedro 



Entanglement in phase space 



56 



Rios also played important roles in the development of the conceptual view that I have 
here attempted to sketch. Maria Carolina Nemes contributed some helpful criticism of 
the manuscript and I also thank the MPIPKS-Dresden for its hospitality during part 
of its development. Partial financial support from Millenium Institute of Quantum 
Information, PROSUL and CNPq is gratefully acknowledged. 



Bibliography 

Schrodinger E 1926 Ann. der Physik 81 109. 
Schrodinger E 1935 Naturwiss. 48 807, 49, 50 844 
Einstein A, Podolsky B and Rosen N (1935) Phys. Rev. 47 777 
Bell J S 1964 Physics 1 195 (see also 0). 

Bell J S Speakable and unspeakable in quantum mechanics (Cambridge University Press, 1993) 
Nielsen M A and Chuang I L Quantum Computation and Quantum Information (Cambridge 

University Press, 2003) 
Bohm D Quantum Theory (Prentice Hall, New York, 1951) 
P. W. Schleich 2001 Quantum Optics in Phase Space(Wiley-VCH, Berlin). 
Bell J S 1981 J. de Physique, 42 C2 41 (see also [5]). 
Feynman R P 1948 Rev. Mod. Phys. 20 367 

Schulman L S 1981 Techniques and applications of path integration (Wiley: New York) 
Rayleigh JWS The Theory of Sound (Dover Publications, 1976) 
Arnold VI 1978 Mathematical Methods of Classical Mechanics (Springer, Berlin) 
Peres A Quantum Theory: Concepts and Methods (Kluwer Academic, Dordrecht, 1993) 
Clauser J F, Holt R A, Home M A and Shimony A 1969 Phys. Rev. Lett. 23 880 
Goldstein H 1980 Classical Mechanics, 2nd edition (Addison- Wesley, Reading, M.A.) 
Voros A 1976 Ann. Inst. Henri Poincare 26 31 

Gutzwiller M C 1990 Chaos in Classical and Quantum Mecahnics (Springer, New York) 
Ozorio de Almeida A M 1988 Hamiltonian Systems: Chaos and Quantization (Cambridge: 

Cambridge University Press) 
Maslov V P and Fedoriuk M V 1981 Semiclassical Approximation in Quantum Mechanics (Reidel, 
Dordrecht) 

Van Vleck JH 1928 Proc. Natl. Acad. Sci. USA 14, 178. 
Ozorio de Almeida A M and Hannay J 1982 Ann. Phys. 138 115 
R. J. Glauber, Phys. Rev. 131(1963) 2766 

J.R. Klauder and B. Skagerstam (1985) Coherent States (World Scientfic, Singapore). 
A. Perelomov (1986) Generalized Coherent States and their Applications (Springer, New York). 
Messiah A 1961 Quantum Mechanics (North Holland: Amsterdam) 

C. Cohen Tannoudji. B. Diu and F. Laeoe (1977) Quantum Mechanics (Wiley, New York). 
Englert B-G, Sterpi N and Walther H 1993 Opt. Commun. 100 526 
L. G. Lutterbach and L. Davidovich,Phys. Rev. Lett. 78, 2547 (1997). 

Bertct P, Auffeves A, Maioli P, Ornaghi S, Meunier T, Brune M, Raimond JM and Haroche S 

2002 Phys. Rev. Lett. 89 200402 
Littlcjohn R G 1995 in Quantum Chaos: Between Order and Disorder, edited by Casati G and 

Chirikov B (Cambridge University Press, Cambridge) 343 
Amiet JP and Huguenin P 1980 Helvetica Physica Acta 53 377 
Coxcter HSM 1969 Introduction to Geometry (Wiley, New York). 
Ozorio de Almeida A M 1998 Phys. Rep. 295, 265 
Grossmann A and Huguenin P 1978 Helvetica Physica Acta 51, 252 
Balazs NL and Jennings BK 1984 Physics Reports 104 347. 
Wigner E P 1932 Phys. Rev. 40 749 



Entanglement in phase space 



57 



[38] Schrodinger E 1930 Proc. Pruss. Acad. Sci. 19 296. 

[39] Deans S R 1983 The Radon Transform and Some of Its Applications (New York: John Wiley & 
Sons) 

[40] Grossmann A 1976 Commun. Math. Phys. 48 191 

[41] Royer A 1977 Phys Rev, A 15 449 

[42] Ozorio de Almeida A M, Vallejos O and Saraceno M 2004 J. Phys. A 38 1473 and 

|quant-ph/ 0410129] 

[43] Zurek W H 2001 Nature 412 712 

[44] Hudson R L 1974 Rep. Math. Phys. 6 249 

[45] Tatarskii V I 1983 Sov. phys. Usp. 26 311. 

[46] Ozorio de Almeida A M and Brodier O 2004 J.Phys. A 37 L245 

[47] Gronewold HJ 1946 Physica 12 405 

[48] Wootters W K 1987 Ann. Phys. NY 176 1 

[49] Miquel C, Paz J P and Saraceno M 2002 Phys. Rev. A 65062309 

[50] Rivas A M F and Ozorio de Almeida A M 1999 Ann. Phys. NY 276 223 

[51] Rivas A, Saraceno M and Ozorio de Almeida A M 2000 Nonlinearity 13 341 

[52] Arguelles A and Dittrich T 2005 Physica A 356 72 

[53] K. Husimi 1940 Proc. Phys. Math. Soc. Jpn. 22, 264. 

[54] K. Takashi 1986 J. Phys. Soc. Jpn. 55, 762. 

[55] Leboeuf P and Voros A 1990 J. Phys A23, 1995 in Quantum Chaos: Between Order and Disorder, 

edited by Casati G and Chirikov B (Cambridge University Press, Cambridge) 507. 

[56] Wehrl A 1978 Rev. Mod. Phys. 50 221 and 1979Rep. Math. Phys. 16 353. 

[57] Bengtsson I and Zyczkowski K 2006 Geometry of Quantum States (Cambridge University Press, 
Cambridge) 

[58] Mintert F. and Ziczkovsky K 2004 Phys. Rev. A 69 022317. 

[59] Mintert F., Carvalho A R R, Kus M and Buchleitner A 2005 Phys. Reports 415 207. 

[60] Bennet H B, Bernstein H J, Popescu S and Schumacher B 1996 Phys. Rev. A 53 2046. 

[61] Glauber R J 1963 Phys. Rev. Lett. 10 84 

[62] Sudarshan E C G Phys. Rev. Lett. 10 277 

[63] Cirone M A 2005 Phys. Letters A 339 269. 

[64] Bell J S 1986 in New Techniques and Ideas in Quantum Measurement Theory (New York Academy 

of Sciences) (see also [5]). 

[65] Banaszek K and Wodkiewicz K 1998 Phys. Rev. A 58 4345. 

[66] Banaszek K and Wodkiewicz K 1999 Phys. Rev. Letters 82 2009. 

[67] Angelo R M, Vitiello, S A, de Aguiar M A M and Furuya K 2004 Physica A 338 458 

[68] Moyal J E 1949 Proc. Camb. Phil. Soc. Math. Phys. Sci. 45 99 

[69] Ozorio de Almeida A M 1984 Rev. Bras. Fis. 14 62 

[70] Tegmark M and Shapiro H S 1994 Phys Rev E 50 2538 

[71] Levy P Theorie de l'Adition des Variables Aleatoires (Gauthier-Villars, Paris, 1954) 

[72] Ozorio de Almeida A M, Magalhaes A and Nemes M C 2006 (in preparation) 

[73] Brodier O and Ozorio de Almeida A M 2004 Phys. Rev. E 69 016204 

[74] Giulini D, Joos E, Kiefer C, Kupsch J, Stamatescu 1-0 and Zeh H D 1996 Decoherence and the 

Appearance of a Classical World in Quantum Theory (Springer, Berlin) 

[75] Diosy L and Kiefer C 2002 J. Phys. A 35 2675 

[76] Agarwal GS 1971 Phys. Rev. A 4 739. 

[77] Feynman RP and Vernon FL 1963 Ann. Phys. (NY) 24 118. 

[78] Caldeira AO and Leggett AJ 1983 Ann. Phys. (NY) 149 374; E 1984 153 445. 

[79] Kenfack A, Rost J M and Ozorio de Almeida A M 2004 J. Phys. B 37 1645 

[80] Carvalho A, Kenfack A, Ozorio de Almeida A M, Rost J M and Toscano F 2006 (in preparation) 

[81] F.A. Berezin and M.A. Shubin (1972) in: Colloquia Mathematica Societatis Janos Bolyiai (North- 
Holland, Amsterdam) 21. 



Entanglement in phase space 



58 



[82] Berry M V 1977 Phil. Trans. Roy. Soc A 287 237-71 

[83] Berry M V 1991 Chaos and Quantum Physics; Les Houches LII (North Holland, ed. Giannoni 

M-J, Voros A and Zinn-Justin J) 251. 
[84] Ozorio de Almeida A M 1983 Ann. Phys. NY 145 100. 
[85] Ozorio de Almeida A M 1992 Proc. R. Soc. Lond. A 439 139 
[86] Rios PPM and Ozorio de Almeida A M 2002 J. Phys. A 35 2609. 
[87] Osborn T A and Kondratieva MF 2002 J. Phys. A 35 5279 
[88] Ozorio de Almeida AM and Brodier O 2006 Ann. Phys. N.Y. 
[89] Dittrich T, Viviescas C and Sandoval L 2005 (private communication) 
[90] Ozorio de Almeida A M 2003 J. Phys. A 36 67 
[91] Ozorio de Almeida A M 1982 Physica A 110 501 

[92] A. I. Shnirelman, Uspehi. Mat. Nauk. 29, 181 (1974); Y. Colin de Verdiere, Comm. Math. Phys. 
102, 497 (1985); S. Zelditch, Duke Math. J. 55, 919 (1987). 



